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SINKARY 


An existing computer program, used foe predicting the natural frequencies 
end mode shapes of helicopter rotor bladaa, has been refined to i approve program 
accuracy end versatility, the program is based on the HoUer-Hyklestad approach 
adapted for rotating beams* Coupled vertical (out -of -plane) , horizontal (in- 
plane) , and torsional mode characteristics can be determined for a variety of hub 
and blade configurations of practical interest. The resulting program is docu- 
mented by presenting the recursion equations and techniques for determining 
natural frequencies and mode shapes, input data requirements and descriptions 
of various program outputs* The accuracy of the program is demonstrated by com- 
paring computed results with exact solutions to classical problems and experi- 
mental data. 


XNTR0WCT10N 

The calculation of undamped natural frequencies and mode shapes or rotor 
blades is an important first step in analysing the dynamic behavior of helicopter 
rotor systems. Air existing analysis and computer program for determining blade 
movies and frequencies# which has been used by both government and industry for 
a number of years, is described in reference 1. This program also provides data 
on punched cards fox input to the flight simulation program desci rbed rn refer- 
ences 2 and 3* It is based on the Holzei -Kykleatad approach adapted foi rotat- 
ing beams and represents the helicopter blade by a lumped-mass system which 
includes effects of rotary inertia# blade geometric pitch, and inertial and 
elastic coupling among the five degrees of freedom at each mass station (radial 
motion is not included). Shear de tor mat ion, aerodynamic effects, and the built- 
in coning (precone) are not represented, oiaciete spanwise variations of blade 
properties and vai u*us hub configurations are Kwxieltst in the program* A unique 
feature of the program of reference 1 is the ability tv' include the effects of 
support system impedance on the blade natural frequencies and rkxW shapes. The 
capability to represent both a variety ot hub coat igur at ions and support system 
etteeta are teatures not generally tound rn othei computer programs ot this 
type. 


Because ot the program 1 * uniqueness, the present authors have made consul 
erable use ot it. However , experience has identified a number of eiiots and 
1 imrtat ioiw» in the computer program* In anticipation ot future need tor this 
modal analysis pr vxji am and possible tuithei expansion ot its capabilities* the 
analysis has been improved and all known errors rn the program have been col- 
lected. Muvu changes tv' the computat tonal pivveduie involved a different iteia- 
t ion scheme to determine natural t requeue res and a different technique ter cvxm 
put ing mode shapes. Other changes to the program were made tv' reduce the side 
ot the source program and coie requirements needed tor execut ion, increase run 
time efficiency, and expand pivxjiam capabilities b^ increasing the allowed num- 
bers of mass stat ions and rotor rotational speeds. The recuisU'n equat ions. 



relating th* state vectors at adjacent stations, were changed in the new program 
version and are substantiated by the work described in reference 4. 

The resulting analysis is extensively different fro* the original version 
and has significant improvements in program site, efficiency, and versatility. 
The purpose of this report is to discuss the equations and techniques used, 
describe the final program, and present a program source listing. The accuracy 
of the program is demonstrated by correlating predicted tesults with exact solu- 
tions for classical problems and with experimental data. 


SYMBOLS 

Values in text are given in U. S. Customary Units and SI units. The calcu- 
lations were made in U. S. Customary Units. The computer program values and 
related material are given in U. S. Customary Units. Conversion factors for 
these units in the computer output are given in appendix A. 

a,b,c distances along A-, B-, and C-axes 

Bj blade section constant, in** (cm**) 

D horizontal force perpendicular to s-axis, lb (N) 

E Young's mxiulus of elasticity, lb/in* (N/»*) 

El blade structural stiffness in bending, lb- in 2 (N-m 2 ) 

F centrifugal force, lb (N) 

in-plane force due to centrifugal force acting through a center of 
gravity horizontally offset from pitch axis, lb (N) 

F x in-plane moment due to centrifugal force acting thtough a center of 

gravity hoiizontally offset from pitch axis, in-lb (N-m) 

Fy out-of-plane moment due to centrifugal force acting thtough a center 

of gravity vertically offset from pitch axis, in-lb (N-ra) 

GJ blade structural stiffness in torsion, lb-in* (N~m*) 

3 acceleration due to gravity, in, sec* sec*) 

1 segment second mass moment ot inertia relative to pitch axis with 

rotary inertia effects included, in-lb-sec- (ro-N-sec*) 

I* cross-sect ional second mass mtvment ot inertia relative to center of 

gravity, lb-sec* (N-see*) 

JHUB number ot nonfeat her ing hub segments 

K 0 blade pitch control system stiffness, in-lb tad {N-m rad) 
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hub in-plaM translational at If (mas par blada, Ib/in (N/m) 

hub out-of-pi am translational stiffnass pat binds, lb/in (N/a) 

rotor shaft torsional stiffMss, in-lb/rad (N-m/rad) 

blada flapping or rotor tar taring spring rata, in-lb/rad (N-A/rad) 

blada lagging spring rata, in-lb/rad (N-A/rad) 

ssgnant naaa radius of gyration, in. (a) 

vortical tores parpandicular to l-axis, lb (N) 

out— of —pi aM banding Aonant, in-lb (N-a) 

hub in -pi ana translational inartia par blade, lb-sac 2 /in. (N-sec 2 /*) 

hub out -of -pi am translational inartia par blada, lb-sec 2 /in. (N-sec 2 /*) 

segment nass, lb-#ac 2 /in. (N-sac 2 /*) 

number of saro frequency modes 

in-plane bending moment, in-lb (N-m) 

rotor rndius, in. (m) 

canter -of -gravity offset from blade pitch axis, in. (m) 
shear canter offset from blade pitch axis, in. (*) 
torsional moment, in-lb (N-m) 

torsional Inertia term representing centrifugal stiffening, in-lb-sec 2 
(N-m-sac 2 ) 


W blade tip weight, lb (N) 

w blade segment weight per unit length, lb/in. (N/m) 

X,Y,5,A,B,C coordinate axes (see fig. 1) 
x,y,s distances along X-, Y-, and 8-axes 
i blade segment length, in. (m) 

s radical distance to end of blade segment, in. (m) 

elastic bending rotat ion in vertical (YS) plane, rad 
A determinant of boundary condition coefficient matiix 
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^ elastic deflection in x-direction (see tig, 1), in* (m) 

tiy elastic deflection in y-directton (see fi 9 . in* (*) 

0 geometric pitch angle between blade principal structural axis and 

hor iiontal (XI) plane at a segaent and, deg 

8 geometric pitch angle between blade principal structural axis and 

hociaontal (XI) plane at a segaent aid-point, deg 

rats# of change of geoaetric pitch angle with span at a segaent aid- 
point, rad/in. (rad/e) 

& c blade root collective pitch angle, deg 

blade linear twist, deg 

X auxiliary function associated with A 

♦ elastic torsional rotation in vertical (XY) plane, rad 

4* elastic bending rotation in horisontal (XX) plane, rad 

ft rotor angular velocity, rad/sec 

<*> blade natural frequency, rad/sec (cpm in tables) 

Subscr iptst 

a,b,c rotating coordinate system fixed to local structural axes at each 
station after deformation 


i, j,l.*A 

or N ♦ 1j 
x,y,« 


blade station number 

rotating coordinate (reference) system fixed to hub 


METHOD OF ANALYSIS 

The method of analysis considers the rotor blade as a discret i«ed system 
of finite elements with each element representing a segment of the blade. The 
user -spec i t ied input data pertains to discretixed structural properties which 
are average values over each segment. The data are referenced to the local seg- 
ment axis system and to either the segment center of gravity (Inertial proper- 
ties) or pitch axis (elastic properties). Since the pitch axis is the reference 
axis in the analysis, the inertial properties are transferred to the pitch axis. 
All properties are then rotated into a hub-fixed axis system oriented so that 
the vertical axis is parallel with the axis of rotation. Details of the data 
conversion are presented in appendix B. Figure 1 illustrates the local segment- 
fixed axis system (A,B,C) and the hub-fixed axis system (X,Y,2) with sign conven- 
tions. The jth segment is shown in figure t having a width ij and an inboard 
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station Sj •• measured fro* the center of rotation. Figure 2 illustrates the 
shear center and center-of-gravity offsets free the pitch axis in the (h,B,C) 
and (X,Y,S) systaais and the sign conventions. The positive directions for r*,, 
r c , and r x are opposite to the B-, C-, and X-axes, respectively. 

The blade segments are converted to a finite-element representation with 
all the inertia lumped at the inboard end and a stiffness element extending over 
the segment length. By using the sign conventions for the displacements, forces, 
and momenta shown in figure 3, using the free- body diagrams of figure 4, assum- 
ing harmonic motion, and requiring force and moment equilibrium and continuity 
of displacements across the lumped masses, the following recursion equations 
result i 
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( 10 ) 


T j " *j*l ♦ “ «y,j+l> ♦ (<•*!««, j ♦ 

♦ [<** ♦ #><«,) ♦ [«*(««» 

whwra th* coefficients in these equations are defined in Appendix B. 

Equation* (1) to (10) are in general agreement with the equation* of ref- 
erence (4), which presents a detailed development of these recursion formulas. 
There are two significant differences* however* between the equations in this 
report and the final equations of reference 4. The equations given in this 
report are based on the assumption that cross-sectional principal structural 
stiffness and mass axes are parallel. Therefore, the terms in the equations of 
reference 4 that arise by assuming they are not parallel are omitted in equa- 
tions (1) to (10). Secondly, the term Fjk| j in equation (1) is not retained 
in reference 4. 

Equations (1) to (10) are used to calculate deflections, slopes, forces, 
and moments at the inboard end of each segment because of the conditions at the 
outboard end of that segment. The calculation process begins by assuming a 
unity value for one of the tip displacements and sero for all others. The cor- 
responding tip forces and moments are not all sere because a lumped mass is 
allowed at that station (tip mass). Thus, for a unity value of the tip out-of- 
plane slope, the forces and moments are calculated from equations (6) to (10) 
and yield 
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where t>, M, and Q represent the out-of-plane slope and moment and in-plane 
moment, respectively, with sign conventions shown in figure 1. Ml other forces 
and moments at the tip are tero. By successively applying equations (1) to (5) 
to determine deflections and slopes at each mass station and equations (f>) to 
(10) to calculate the forces and moments acting at the outboard end of the next 
inboard elastic element, the conditions all along the blade are determined for 
specified values of collective pitch angle and blade twist, rotor rotational 
speed, and the assumed trequency of vibration. This procedure is repeated for 
the same vibratory frequency by individually assuming unity values for the other 
tip displacements (4 N+ i, $ X| h+ 1, $y,na1>. 

To analyte a particular rotor blade, it is necessary to specify the geo- 
metric and structural properties relating the manner in which the blade is 
mounted to the hub* The hub configuration and support system impedance charac- 
teristics are expressed as five boundary condition equations in terms of the 
displacement variables, forces, and moments at a particular station, IV satist^ 
the prescribed boundary condition equation with the conditions calculated along 
the blade by equations (1) to (10), the unknown u' must take on only select val 
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urns. Thus, the eigenvalue problem Is formed where the vibratory frequency is 
systematically varied, and the conditions along the blade are calculated for 
each unity tip displacement. These data are then used to see whether the bound- 
ary conditions are satisfied. 

In the analysis three basic sets of boundary condition equations are used 
and the resulting nodal characteristics are designated as pertaining to collec- 
tive, cyclic, or scissors type nodes, depending upon which set of boundary condi- 
tions is inposed. The forn of the boundary conditions for each node type is 
susnar ited in table 1 and the different spring and mass terms are shown in fig- 
ure 5. Reference S presents a discussion of these mode types and how they are 
combined to describe rotor mode characteristics. 

The collective mode is characterised by symmetric vertical or out-of-plane 
and antisymmetric horlsontal or in-plane deflections of opposing pairs of blades 
on the rotor. The appropriate boundary conditions for the in-plane direction 
are elements 2 and 4 of the first column of table I and pertain to a spring- 
restrained K»c pinned joint at the center of rotation. The first and third 
elements give the conditions for a clamped joint attached to a movable hub, 
where the hub impedance is represented by a sinqle-degree-of-freedcm mass-spring 
system. These characteristics describe the boundary conditions for the out-of- 
plane direction. The first four conditions of the first column are applied at 
the center-line station. The torsional equation relates the twist and torque at 
the pitch horn radial attachment point. The term represents the effective 

spring rate of the control system. 

The cyclic modes have symmetric in-plane and antisymmetric out-of-plane 
deflection shapes about the center of rotation. The boundary conditions for 
the cyclic modes are given in the second column of table I. The first and third 
elements of that column describe the flapping-spring-restrained pinned con- 

ditions in the out-of-plane direction. Elements 2 and 4 are for the in-plane 
dnection where the representation is a clamped joint fixed to a flexible hub 
having lumped stiffness and mass char act • r i sties described by Kjp and Mjp, 
respectively. The torsional boundary condition is the same for the cyclic and 
collective modes. 

For the scissors modes, the in-plane and out-of-plane boundary condi- 
tions at the center line represent clamping of the blades to an immovable hub 
(\,1 • | • • <l’i • 0). For a rotor having three or more blades, the tor- 

sional condition is the same as for the collective and cyclic modes. For two- 
bladed rotors, the torsional deflection is zero (clamped) at the pitch horn 
attachment radial station. These conditions are shown in the third column of 
table I. For the in-plane and out-of-plane directions, an alternate form of 
the boundary conditions is used whenever it is desirable to represent offset 
flapping and lagging hinges as in the ease of an articulated rotor . (See the 
fourth column of table I.) If the offset of the flapping hinge is not zero, 
the *ero slope condition is replaced by an equation relating the out-ot-plane 
moment and slope at the hinge station using a flapping spring term . Simi- 
larly, for an offset lagging hinge, the slope condition is changed to a model 
of a spr ing-restrained Ity pinned joint at the lagging hinge radial station. 
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Tha valuta of deflection, slope. Moment, and shear calculated at the sta- 
tions for which the boundary conditions apply and due to a value of unity for 
one coordinate at the outboard tip can be substituted into t’ * left-hand side 
of the boundary condition equations for a particular mode type. These equations 
are then uaed to form one column of the boundary condition coefficient matrix. 

By repeating the substitution of conditions associated with unity values '"'e 
other tip coordinates, the complete coefficient Matrix is generated, t. M s,i . 
terms may be thought of as partial derivatives of each boundary condi ♦ .1 equa- 

tion with respect to an individual tip displacement. As an example, to; a Mad* 
having a pitch horn offset the coefficient matrix [c 3 for the scissors mode 
would be (see table X) 1 
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36 y»i 
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3 fi y,i 

3 «y,l 


3e tip 

3,J, tip 

3 ^x,tip 
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3 ^x,l 

3f x,l 
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3 4*tip 
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3tS tip 
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3ik<i x, tip 

34'tip 

^y.tip 


3(T - K c $) 
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the subscript 
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to the station 

at the center 

of rotation* 

The 


form of the boundary condition coefficient matrix is dependent on the mode type 
to be calculated. Because the boundary condition equations are homogeneous, the 
following matrix equation can be formed to determine the tip deflections for 
each mode: 
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Fc»r nontrivial values of tip deflection, the determinant of the coefficient 
matrix must vanish and yield a polynomial in terms of the squares of all natural 
frequencies. In general, the variation of the determinant of the matrix with 
frequency is such that only iteration techniques, where the frequency is incre- 
mented to find a crossover point, are reliable. Such frequency stepping tech- 
niques, however, are inefficient. As shown in reference h, an auxiliary func- 
tion can be generated which obviates the need for frequency stepping. The 
auxiliary function has the desirable feature of monotonically decreasing with 
increasing frequency for frequencies less than the natural frequency being 
sought. (See fig. $.) The behavior of the auxiliary function is du*» to the 
removal of all toots of the determinant below that of interest and allows 
use of extrapolation techniques for convergence to each natural frequent e 
auxiliary function is defined by the relation 


\i («') 
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to* - .OJ*! fo* - .0^*) . . . 
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where the ith natural frequency is desired. At the ith natural frequency, the 
values of \ and A are both se.o. In the computer program the auxiliary func- 
tion is used to extrapolate to tire natural frequencies. Once the natural fre- 
quency is determined, that toot is also removed t tom the determinant before pro- 
ceeding with the calculations for the frequency of the next mode. In practice, 
the P »ero frequency roots, which correspond to the rigid- tody modes, must al so 
be removed before this technique will wotfc p'ropvrly (a fact not mentioned in 
ref. hi. This removal may be aocomp.l ished by noting the slope ot the auxiliary 
function feu values of frequency near sere. If the magnitude of V increases 
for increasing frequency, then \ is divided by successively higher powers ot 
,e* at the trial point until convergence is assured. The .o* p factor is main- 
tained in the denominator of the auxiliary function leg. (HI) for all successive 
mode calculations. 


Subst i tut ion of the calculated natural quency into the tvundaty condi- 
tion coefficient matrix leads to five hcmog< •. o.rs equations in terms ot the 
fire unknown displacements .•,( the outtvard trp (eq. (l.‘l. To solve tor the 
relative magnitudes ot the tip' deflections, an inverse iteration technique 
(ref. ?l is used. The inverse of the ooet t tor ent matrix is used to ptemultvply 
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•n initial trial vector. The resulting vector is used as the new trial vector 
and is again premultiplied by the inverse matrix. Pour iterations are perforated, 
each resulting vector being normalised by the largest element. The final vector 
corresponds to the relative tip deflections for that node. By using these val- 
ues and the spanwise distributions of deflections, slopes, shears, and moments 
for each unit tip deflection, the node shapes and associated shear and moment 
distributions are calculated. The displacements are left in the X,Y,Z-axis 
system. The shears and moments are resolved back into the loca segment axis 
systems, however, by using the twist and collective pitch angles. 


DESCRIPTION OP COMPUTER P:<OGRAN 
General Description of Program 

The analytical methods have been implemented as a digital computer pro- 
gram coded in FORTRAN IV language. The program is run on the Control Data 
Gorp. CYBER 175 computer system at the Langley Research Center, using an FIN 
compiler and NOS 1.2 operating system. By comparing the new program with the 
one described in reference 1, the lengths of the source decks are 1222 and 2398 
cards, compilation times are 3.6 and 66 CPU (Central processing unit) seconds, 
and computer core requirements in octal are 54000 and 112000, respectively. A 
sample case with a 20-regment representation was run or. both programs. The new 
version required 25 CPU seconds for execution and the oid version required 27 
CPU seconds. The new program can calculate natural frequencies and mode shapes 
for up to 30 blade segments (user specified in the new program), and up to 10 
rotor speeds. All errors in the original program known to the authors have 
been corrected in this program. 

The computer program consists of a main program, called BLDANL, and eight 
subroutines (ANPLTD, CARDS, COST, INPT, PLOUT, START, SUMMY, and MODS HAP) . A 
listing of the computer program is given in appendix C. In uddition, a library 
subroutine, MATINV, is also required and is documented in appendix D. The pro- 
gram also uses a CalCamp plotting package. A flow chart of the program is pre- 
sented in figure 7. The various parts of the program are described in subse- 
quent paragraphs. 


Program BLDANL 

The main program controls the computation of the blade natural frequencies 
for all combinations of collective pitch, rotor speed, and mode type. BLDANL 
initially calls subroutine INPT to read the input data for each case and then 
convert the input into the desired lumped-mass representation. The elastic and 
inertial properties for each station are rotated from the local A,B,C-axis 
system to the hor i son tal- vertical X,Y, 2-axis system. The program then computes 
the forces due to rotation at each station for each value of rotor speed. 

The calculated characteristics in the X,Y,Z-?xis system are input to 
subroutine COBP which then returns the determinant of the boundary coefficient 
matrix (eq. (12)) for each desired estimate of the natural frequency. In prac- 
tice, the values of the determinant are used to extrapolate toward the point 
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«h«* the determinant ha* a itto value. Once the program has converged on the 
first natural frequency, the process la continued for successively higher nodes 
until the specified upper frequency is exceeded. This procedure is followed 
far each appropriate node types collective, cyclic, or scissors. 

After all desired natural frequences have been determined, the associated 
node shapes nay be calculated. A call to subroutine AMPLTD causes the node 
shapes to be calculated, printed out, plotted, or punched depending on the input 
options selected. Subroutine SIMOT is then called to print out a suanary of all 
calculated natural frequencies. Also, if requested, subroutine FLOUT is called 
to plot the natural frequencies as a function of rotor rotational speed for each 
collective pitch angle and node type. The prog ran then proceeds to the next 
input case. 


Subroutine AMPLTD 

Subroutine AMPLTD is called by BLDANL to compute the node shapes for the 
previously determined natural frequencies. Subroutine OUST is first called to 
calculate the matrix of the root-boundary conditions for individual unit tip 
displacements at the calculated natural frequency. The matrix is inverted by 
using the library subroutine MATINV and an inverse matrix iteration scheme is 
used to calculate the values of the deflections and slopes at the tip. From 
the transformation terns computed in subroutine OOEF, the deflections, slopes, 
moments, and shears at each station are calculated and listed. The part of the 
generalized mass associated with each deflection or rotation cor' went is calcu- 
lated and used to determine the principal deflection direction (i>at having the 
largest generalized mass): vertical, horizontal, or torsion. The uumber of 

node points of the mode shape in the principal direction is then determined. 

If requested, the mode shapes may be plotted or punched out through calls 
to MDDSHAP or CARDS. The generalized mass components are also printed. Finally, 
a check of convergence for both natural frequency and mode shape is made tv- 
matrix multiplication of the boundary condition matrix and the matrix of tip 
deflections and slopes. The retulting matrix, which should be a null matrix, 
is printed out for each listed mode shape. 


Subroutine CARDS 

The subroutine CARDS is the only subroutine which produces punched output 
for use as input data in the flight simulation program of references 2 and 3. 

It is divided into three parts - the first punches the input inertia data, the 
second punches the components of the mode shape, and the third punches the 
cyclic detuning cards which ore used to specify variations of natural frequency 
with collective pitch and rotor speed. The first part, called by subroutine 
START, takes blade mass and inertia data which are divided into an arbitrary 
number of segments and recasts it into 20 equal length segments. The second 
part, called by the subroutine AMPLTD, takes the mode shape data and recasts 
it into 20 equal segments as well. The first six natural frequencies for each 
combination of rotor collective pitch, rotor rotational speed, and mode type 
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ars stored. These frequencies «r* used in th* third part to compute informa- 
tion for th# cyclic do tuning cards. 


Subrout in* COW 

the aubroutina COST is used to for* th# root boundary condition coefficient 
Matrix tor a apaclf i*d frequency, Th* fiva generalised coordinates at th* out- 
board tip - vortical slope, horizontal slop*, horisontal deflect ion, torsional 
rotation, and vortical deflection - ar* individually sat to a valua of unity 
whereas th* restsining tip daflactions and alopaa ar* sat to raro. by using th* 
alastic and inartial properties or tha X,Y,i-eoordinat* system and th* speci- 
fiad frequency, tha deflection, rotation, shear, and moment at each station along 
th* blade ar* calculated for each unit tip displacement. Fur the pai t icular 
node type and rote* configuration the combination of five root condition* to be 
zeroed are selected. The previously calculated values of the condit ions on the 
blade at the boundary station# are inserted into the appropriate row and column 
of the coefficient matrix. When the boundary coefficient matrix is eoatpleted, 
aubroutine HAT1NV is called to compute the determinant of the matrix or its 
invar a*. 


Subroutine INPT 

The subroutine INPT is called by BLDANL to read the input data from punched 
cards. INPT first reads the program option card and then either a full data 
dec* or changes to the previous case. For multiple-case runs, a namelist format 
may be used for all cases except the first. If desired, subrout Ine CARDS may be 
called to output the blade inertia properties on punched cards. Subroutine 
START is also called trim INPT. 


Subrout ine WttiT 

The aubroutine PLOW produces plots of the variation or natural frequencies 
with rotor speed, fine figure is produced tot each type of mode. For each fig- 
ure, data for all eombinat ions of collective pitch angle ami rotiH rotational 
speed are plotted. The maximum Inertial pl.w Associated with each natural fre- 
quency ia distinguished by different plotting symbol*. Subroutine FLOUT is 
called by program BLDANl.. 


Subroutine START 

This subroutine is called by INPT. The location* of the pitch horn. Hap- 
ping, ami lagging hinge offsets are determined in tetma ol segment number and 
location within the appt opt late segment. The centrifugal force act mg at each 
station is calculated along with blade mass and second nw'ment ot ine» t ia 

alxxit the flapping hinge or center line. Km linear twtst dist t ibut ions, the 
twist at each station is also calculated in subtout ine START. The input inettra 
data which pet tain to individu.il segments ate averaged to generate inertia prop- 
erties at each station amt ptoducta and or oss> pr oduct s ot met t ia ate calculated 



ter **ch station m mU. tbm Input data ut printed out (or identification 
purpoMt. 


Subroutine SUNK 

Subroutine SIAM? is called by the Min progras BLDANL to list a suamary of 
the natural frequencies for each combination of rotor collective pitch, rota- 
tional speed, and sods type. Listed along with the natural frequencies are the 
maximum inertial (principal deflection) plane pertaining to each frequency and 
the number of node points of the component of the mode shape identified by the 
inertia plane. 


Subroutine NOOSHAP 

Subroutine NOOSHAP is called by subroutine AHPLTD to plot out the mode 
shape. The hoc ixontal and vertical deflection and the torsional rotation are 
plotted as a function of the blade radial station nondimensionalised by the 
radius. Sad) component is identified by a distinct symbol. For convenience, 
the associated natural frequency is also given in the plot. 


Description of Program Input Requirements 

The program reads the input from file 5. The format for the input data 
deck is specified in table I1 and a sample deck listing is provided in table 111. 
The program input and output oata are expressed in U.S. Customary Units. Conver- 
sion factors for obtaining SI Units are given in appendix A. table II the 
program names are given with the symbols used in this report shown parentheti- 
cally, where appropriate. A data deck consists of card types 1 to 18. Multiple 
cases can be run by placing the data decks one behind the other with no separa- 
tors. For multiple-case runs, the namelist forMt option can be used for all 
cases except the first. The various inputs are described in the following 
paragraphs. 

In table II, the first card apecifies the program options via descriptive 
names which may be input in any order. The allowed names are DECK, NAMElist, 
PUNCh, SHAPe, NODES, ALLModea, PLOT, TQRSion, and NLTNist. Only the first four 
characters o' each name are required. DECK instructs the program to read a new 
case input deck. The NANElist option allows the user to read in only the changes 
to the previous case by using a namelist format. PUNCh is used to generate 
punched cards suitable for input into the computer program of references 2 and 
3. The SHAPe option causes the program to generate CalComp plots of the blade 
mode shapes for the case with the reference values of rotor speed and collective 
pitch angle. The NODES input results in printouts of the mode shape data pro- 
vided by SHAPe. ALLModes is used when mode shape printouts for all combinations 
of rotor speed and collective pitch angle are desired. The PLOT option gener- 
ates CalComp plots of the variation of blade natural frequencies as a function 
of rotor speed for each collective pitch angle. TORSion causes the torsional 
degree of freedom to be included in the analysis. The NLTOist option instructs 
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the proftHM tg i««d in * nonlinear tmi»t distribution uaing card 11 (card 1 ) i* 
emitted It NtdMiat U net u*ed). 


Card type 3 i* treed to th# restraint condition* imposed on th* 

bled* notion* by th* pylon and collective pitch control system*. Th* pylon 
restraint U provided by lumped ****** and spring* in th* out-of -plane and in- 
plan* dir action* two*, HIT* HOP, and KIP) and a m**t toraional wind-up spring 
KT, Tha control *y*ta« *titfna«* i» represented by KC. Tha tarn KWCTA nay be 
ueed to proa id* a rotor teetering apt in* or blada flapping spring in th* caae 
ot an articulated rotor* Tha lagging spring rata ia input by KPPl. 

Tha blada nay ba dividad up to a maximum ot 3b segment*, tha inboard d»fU» 
segment* bain* nonfeather in*, it all segments ar« ot egual length, a nonaero 
value of AKMK ia input and card typ* 7 ia omitted. Nonuniform segment length* 
ar# input by tatting atftJdt to *a> and apacitying tha radial station of tha 
outboard and of each aagnant on card ?. Th* rotor geometry ia apacifiad by tha 
blada WIST or TWO* number ot PUtbJS, and tha radial location* of tha lagging 
tCNOPP) and flapping (INOPf) hinges, and pitch horn (PKOPP). to estate spe- 
cif ie hinge*, a taro input valua ia uaed, The tarn ia uaad only to aup- 

praaa tha calculation ot collactiva and cyclic node data when tha input valua 
i* greater than 2.0, 

Tha program datarninaa natural treguencies bale* tha apacifiad upper Unit* 
PIAST, for all combination* ot rotor collective pitch angle ikcou.) and rotor 
apaad tNftPM), 1%' to 16 rotor speeds U*m) and i collective pitch angle* 

(IftCQh) may be invaatigatad tor any one case, However, it tha natch option ia 
ueed, three collaotiva pitch angle* ami three rotor apaeda nuat ba input, aa 
required by tha program of rererence* 2 and 3, tha second value* of each param- 
eter being tha average ot tha fit at ami third value*. The SHAPe ami MOPS* 
option* provide output data pertaining to tha next to the higheat valua* ot 
both rotor speed ami collective pitch angle if three or mv« value* are ape *- 
tied, for one .u two value* ot rotor speed ami or collective pitch, the refer- 
ence value* are the first ot second, r eapeot ively, 

Kach segment is assumed to have oonstant properties over its length. Card 
type* f) to 10 ami 12 to 1* are used tv' input average segment charaotet iat ica, 
eight values pet card, card 10 the tip weight ia input after the N values 
vvf segment weight pet unit length itmn.l. The mass moments ot inertia IKYKB 
and KtKC) are resx'lvevl about rhe center of gravity. ,- Se locat rv'ns ot the shear 
.'enter ami .'enter of gravity relative to the pitch ax*,* are input based on the 
sign conventions ot figure 2, Note that the positive directions tor Kb and Hr' 
{center -ot-gr av ity offsets rn b- ami c-driect tons) are opposite to the b~ ami 
C-axi* convent ion. Tire N ♦ 1 value ot Kb and IK' pertain to rhe eccentricity ot 
the tip weight relative to the pitch axis. 


Pescrrpt ion ot t'r.Mraw Output l\*ta 

The tntormatton output by the pr.xjram rm'ludes printer listings ami, 
optionally, punch#-.. cards .u Cal.'.xxp plots depending .ur the program .vntrcl 
»yt ion* selected by the user .m\ card tyr>e 1. iSee discuss r.\n .u\ input data 
r egu tr ement s, 1 The pmvgtam generates a listing ot the ident it toat ton, geo- 
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»«tric, and structural data that were input and a smeary table of all calcu- 
lated natural frequencies. Samples of these outputs, obtained fraa the input 
data listed in table III, are presented in tables IV and V. 

The input data list includes geometric and structural parameters in the 
same units as originally read in, as well as segment lengths and values of blade 
twist and centrifugal force at each station. The total blade mass and flapping 
inertia are calculated by considering only the blade mass outboard of any flap- 
ping hinge and are listed. The hub and control system impedance parameters are 
also reproduced on the listing. The calculated natural frequencies are summa- 
rised for each appropriate mode type (that is, collective, cyclic, and/or scis- 
sors) and at each combination of root collective pitch angle and rotor speed. 

The valuis of natural frequency are nondimensionalised by rotor speed prior to 
being output or they are presented in units of cycles per minute for nonrotat- 
ing cases. For each frequency listed, the plane containing the largest contri- 
bution to the total generalised mass of that mode is identified and the number 
of node points in the mode shape component associated with that plane are 
listed. 

If either the MODES or ALlModes options is selected, spanwise distributions 
of the mode shapes, moments, and shear forces at each station are printed out as 
shown in table VI. When the ALLModes option is exercised, these data are listed 
for all combinations of appropriate mode type, root collective pitch angle, and 
rotor speed. For the NODES option, data are printed far each mode type, but for 
only one reference combination of collective pitch and rotor speed. The refer- 
ence collective pitch angle is the last value input if only one or two input 
values a e specified. If three or more collective pitch angles are requested, 
the reference value is the next to last value input. The reference rotor speed 
is determined in a similar fashion. The values in the mode shape columns corre- 
spond to deflections in the vertical and horizontal planes. However, the asso- 
ciated forces and moments are calculated in the local beam and chord axis system 
and reflect blade built-in twist and collective pitch angles. When TORSion is 
input on card type 1, the torsional rotations and moments for each mode ai e also 
listed. The displacements, moments, and shears are normalized so that the maxi- 
mum defJ ccion along the blade in either the vertical or horizontal plane is 
1 in. (2.54 cm) or the maximum torsional rotation is 1 radian. After the span- 
wise mode shape listings, the parts of the total generalized mass attributed to 
each deflection plane, vertical, horizontal, or torsion, are printed. Finally, 
a check is made to see whether the boundary conditions are adequately satisfied 
and the calculated boundary condition values are printed. (See previous discus- 
sion concerned with imposing boundary conditions.) Both the generalized masses 
and boundary condition values reflect the normalization process applied to the 
miide shap 's , 

U of the SHAPe option causes CalCcmp plots of the spanwise distribution 
of * * mode shapes to be generated. A sample plot is shown in figure 8. One 
plot is made for each calculated frequency for all the mode types but only for 
the one combination of reference collective pitch and rotor speed. The differ- 
ent displacements, vertical, horizontal, or torsion, are plotted individually 
iind reflect the normalization process applied to the listed data (maximum deflec- 
tions of 1 ir« (2.54 cm) or 1 radian). The abscissa represents the spanwise sta- 
tions for :he deflection values normalized by the rotor radius. For identiftca- 
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fcion purpoiM, the mode type and natural frequency for aach plot ara given above 
the figure. 

The PLOT option produces aaparata plots of blade natural frequency as a 
function of rotor spaed for collective* cyclic* and scissors node types. For 
eadi plot* a saaple of which is shown in figure 9* the natural frequencies for 
all combinations of collective pitch angle and rotor speed are shown. The 
plotted data distinguishes between siodas having maximum generalised masses in 
the three component directions (vertical* horisontal* or torsion). For identi- 
fication pv -poses* the case number* title* and other information are given in 
each plot. 

When the FlNCh option is used, punched cards are generated by the program 
compatible with the input requirements of the rotor craft simulation computer 
program documented in references 2 and 3. The format for the punched deck is 
given in table VII and a sample output deck is listed in table VIII. The 
punched output consists of four parts: a problem identification card* cards 

containing the blade inertia characteristics, the mode shapes for the six low- 
est frequencies associated with each mode type and cyclic detuning information. 
When the PCNCh option is used* three values each of collective pitch and rotor 
speed must be input* the second values being the average of the low and high 
input quantities. 

Although the program has the capability of treating up to 30 segments of 
unequal length, the simulation program of references 2 and 3 requires data for 
only 20 serpents of equal length. Therefore* the input rotor inertia character- 
istics are converted by using linear interpolation. The moments of inertia in 
both beam and chord directions are also changed to reflect a transfer from the 
center of gravity to the feathering axis prior to being punched out. Bach set 
of inertia properties is punched out in successive fields on three cards. The 
calculated mode shapes for the reference values of rotor speed and collective 
pitch have also been interpolated to yield values for a system with 20 segments 
of equal length. The units of the mode shape components are changed to feet for 
the vertical and horisontal deflections and degrees for the torsional amplitudes. 
For each of the first collective* cyclic* and scissor modes (up to six each) * 

11 cards are output. The first 10 cards define the three mode shape components 
(6 X 6y* and ♦) for adjacent station pairs starting at the center of rotation. 

The tip displacements are punched in the first three fields of the 11th card 
with the next three fields containing the natural frequency, normalized by the 
rotor speed* mode type designation* and an assumed modal damping ratio of 0.02. 
Identification data are punched in the remaining fields of all mode shapes cards 
for ease of handling. The mode type designations are 1 (collective), -1 (cyclic), 
and 0 (scissors) . 

One cyclic detuning card is generated for each punched mode shape. The 
purpose of the cyclic detuning card is to specify the variation of mode fre- 
quency with rotor speed and collective pitch angle. The first four fields of 
each cyclic detuning card list the normalised frequency of a particular mode 
for conditions of lowest values of rotor speed and pitch angle, lowest rotor 
speed and highest pitch angle, highest rotor speed and low pitch angle and, 
finally* the highest values of pitch angle a^d rotor speed. The reference val- 
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um of rotor spaad and pitch angle are also punched on each cyclic detuning 
card, as well as the case nans, node number, and type Indicator. 


PROGRAM CORRELATION RESULTS 

The accuracy of the analysis is verified through correlation of computed 
results foe selected problems, having either closed-form solutions or accurate 
solutions available in the literature. The correlation is based on comparisons 
of natural frequencies and node shape characteristics. The first case consid- 
ered is that of a nonrotating, uncoupled, uniform beam with pure bending in one 
plane and torsion allowed, (tee-free and clanped-free conditions are examined 
and the lowest five modes in each direction and for each set of boundary con- 
ditions are determined. The second case expands on the first by adding a tip 
mass in a manner which affects the bending vibrations only. The third case is 
for a simply supported, nonrotating beam having uncoupled torsion and one bend- 
ing degree of freedom. The fourth case adds elastic coupling between the bend- 
ing and torsion vibrations for the simply supported beam. The fifth correlation 
case deals with a rotating propeller for which experimental data were available. 
Unfortunately, suitable data pertaining to m helicopter rotor could not be found 
by the authors. The analytical representations for each of these cases are sum- 
marised in tables XX and X. 


Case 1 - Nonrotating Uncoupled Uniform Beam 

The program input data for the nonrotating, uncoupled, uniform beam in 
bending and torsion are presented as case number 1 in table IX. Although the 
exact solution does not include the effects of rotary inertia, the computer pro- 
gram does. Therefore, negative values of the sectional mass m o ment of inertia 
(determined by eq. (Bl)) are input to the program so that the net segment rotary 
inertia effects in the progr.ua are eliminated. Comparisons of computed and 
exact results (refs. 8 and 9) for the first case are tabulated in tables XI to 
XIV. The differences between the computed and exact natural frequencies for 
both sets of boundary conditions are always less than 0.8 and 1.7 percent for 
bending and torsion, respectively. The discrepancies increase with mode number 
and the computed frequencies are always smaller. The trends of bending mode 
shape correlation are similar to those for natural frequency. The computed and 
exact torsion mode shapes match for all nodes and both sets of boundary 
conditions. 


Case 2 - Nonrotating Uncoupled Uniform Beam With Tip Mass 

Table XV illustrates the effect of adding a tip mass on the clamped-free 
bending frequencies. The analytical model is shown in table IX as case num- 
ber two. The torsion frequencies are unchanged from the values indicated in 
table XIV. The calculated bending frequencies vary by up to 1.5 percent rela- 
tive to the values obtained from reference 10 (for the first vertical bending 
mode). The trend of increasing discrepancy with mode order observed for the 
previous examples is not present in this case. In fact, the smallest difference 
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occur* for the third bonding oodo (0.3 poroont). for tho rongo of dot* in 
toblo XV, tho oorrolotlon is still good. 


Cosos 3 ond 4 - Uniform Siaply Supported Boon 

Coupled bonding-torsion vibrations aro analysed in reference 11, where equa- 
tion* are given for calculating exact natural frequencies and node shapes. The 
example problea is that of a unifora, siaply supported (pinned -pinned vertical 
bending and c leaped -claaped torsion) beaa with a horixontal offset between the 
center of gravity and shear center axes to induce coupling. These boundary con- 
ditions cannot be achieved at the outboard tip by the analysis docuaented in 
this report without slight aodifications. The required changes to the ccaputer 
prograa include replaceaent of the cod* generating the unity values of vertical 
one torsional deflections and all teras pertaining to shears and aoaents at the 
outboard tip. Code representing unity values of the vertical shear and torsional 
aoaent are then substituted into the prograa. These changes were aade and the 
analytical representation, given as cases 3 and 4 in table XX, used for correla- 
tion purposes. Tables XVI and XVII present coaparisons for calculated and exact 
results with and without the coupling added. Without the coupling the largest 
discrepancies in frequency are less than 0.02 and 1.7 percent for the highest 
bending and torsion Modes, respectively, thus verifying the representation. For 
the coupled problea (table XVII) the differences between computwd and exact fre- 
quencies are always less than 0.7 percent with the difference generally increas- 
ing with node order. The relative bending- tors ion amplitudes for the node shapes 
show good agreement between computed and exact results. 


Case 5 - Model Propeller 

Experimental vibration data, pertaining to model propeller blades, are 
given in reference 12 and were used in reference 13 for correlation with the 
analysis prasented in that report. The structural characteristics present*-' in 
these two references were used to obtain the representation shown in table X by 
interpolating to intermediate stations. The effective torsional stiffness which 
was used included both the nominal value of GJ and the contribution of another 
term dependent on the square of the twist rate as shown in table X. The second 
term is included in the torsional equation of motion in reference 13 but not in 
the equations given in this report (necessitating its inclusion with the GJ 
values). Me effect of this term is to increase all frequencies alike without 
regard to the values of collective pitch or rotor speed. Its amission in the 
present equations is based on its insignificance for the twist rates used on 
conventional rotor designs. 

Experimental and calculated variations of natural frequency with collective 
pitch angle and rctor rotational speed are illustrated in figure 10 for the 
first and second bending and the first torsion natural frequencies. Excellent 
agreement between experimental and calculated results are achieved for the two 
bonding frequencies. (See figs. 10(a) and 10(b).) A discrepancy exists between 
calculated and measured variations of torsional frequency with collective pitch, 
as shown in figure 10(c). As discussed in reference 13, the angle between local 
structural axes on the blade and the plane of rotation is a function of built-in 
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twist, collective pitch angle, and steady-state twist deflections caused by rota" 
tion. The last item la not included in tha linear analyaia documented in thla 
repeat. However, thia tarn ia not significant for conventional hall copter blade 
designs. 


CONCLUDING REMARKS 

An existing computer program, used for predicting tha natural frequencies 
and node ahapaa of helicopter rotor bladea, haa bean modified to improve program 
efficiency, accuracy, and veraatility. The aquations and techniques used in 
the original version ware extensively changed. Thia program is baaed on the 
Bolter -Myk lea tad approach adapted for rotating beams and represents the rotor 
blade as a aeries of lumped masses connected by stiffness elements. Elastic 
and inertia coupling between vertical, horitontal, and torsional deflections 
are included in the analysis. Spanwiae variations in structural stiffness, 
mass, blade twist, rotary inertia, and center of gravity and shear center offsets 
from the pitch axis may be input. Provisions are also included for representing 
various hub configurations and the effects of pylon impedance on blade modes. 

The resulting program is documented by presenting in thia report the equa- 
tions and techniques for determining natural frequencies and mode shapes, input 
data t'tquir extents, and descriptions of the various program outputs. A source 
listing of the program is also presented. The accuracy of the program is indi- 
cated by comparing computed results with exact solutions and experimental data. 
For uniform, nonrotating beam pro bi«wa, natural frequency predictions are within 
2 percent of exact solutions for the first five b* Jing and torsion modes and 
for a variety of boundary conditions. The predicted node shapes closely match 
exact answers. Experimental data measured on a rotating, highly twisted propel- 
ler were also used for correlation, the calculated bending frequencies showed 
excellent agreement with measured data. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
April 20, 1971 


20 



APPENDIX A 


CONVERSION OP U.S. COSTONARY UNITS TO 81 UNXT8 

The input and output data associated with the computer prog care art 
expressed in U.S. Customary Units. conversion factors for obtaining data 
in SI Units ara givsn in tha following tab la basad on rafaranca 14c 


Physical quantity 

Units 

Convarsion 

factor* 

U.S. Customary 

SI 

Porno, Maight 

lb 

N 

4.4482 

Frequency 

cp» 

Rs 

0.0K6? 

Length . . 

in. 

ft 

0.02540 

Linaar spring 

lb/in. 

N/» 

175.13 

Naas 

lb- sec 2 / in. 

N-sac 2 /» 

175.13 

Monant 

in-lb 

N-m 

0.11298 

Rotational spaad ......... 

i 

rpm 

rad/sac 

0.10472 

Rotational spring ......... 

in-lb/rad 

N-m/rad 

0. 1 1 298 

Second mass moment of inertia . . . 

lb-sec 2 

N-sec 2 

4.4482 

Structural stiffness . • 

lb- in 2 

1 

N-m 2 

0.28698 * 10“ 2 

L 


•Multiply valuas givan in U.S. Customary Unit by convarsion factor to 
obtain equivalent valua in SI Unit. 
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APPENDIX E 


DEVELOPMENT OP TEAMS APPEARING IN THE RECURSION EQUATIONS 

To determine the natural frequencies and and* shapes associated with a 
rotor blade* tha blada it segmented* and average values of tha struct ucal pr op- 
art! aa of aach of tha N segments ara determined. Thaaa proper ti as ara associ- 
ated with tha midpoint of aach segment by tha prograa and tha cross-sectional 
second aaas moments of inartia ara ralatad to tha section aass centroid* Fig- 
ure 1 illustrates tha local A*B*C segaent-fixed axis systea and overall X,Y,S- 
sxis systea with sign conventions* Tha reference axis is tha undetoraed l 
pitch axis of tha blade. The jth segment is shown in figure 1 having a width 
ij and an inboard station Sj as measured from tha center of rotation. 

The aass aoaents of inertia of the segment cross section for the beam and 
chord directions l£jj and 1^ are transferred to the reference (blade 
pitch) axis and tha rotary inertia for the segment is added. The polar moment 
of inartia is also determined by the program. The relationships aret 


» . "i*i/»i 2 2 \ 

ibba • ibb,i*i ♦ — \\T + rb *V 


(i - 1, 2, .... N) <B1) 


I cc,i * J cc,i*i 


"i*i/*i 2 

♦ 1 — + 

g \12 



(i » 1, 2, . . N) (B2) 


*aa»i " t 5 bb,i ♦ J cc,i)*i 



r c,i) 


(i - 1, 2, . . N) (B3) 



I aa,i9 

*i w i 


(i » 1, 2, . . ., N) (B4) 


where w* is the weight per unit length for the ith segment. The terms 
and r c j are the center-of-gravity offsets from the reference axis and the 
sign convention for these terms is given in figure 2(a). 

Because the analysis represents each segment as a weightless stiffness ele- 
ment with the inertial properties lumped at its inboard end, the inertias for 
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M-MSND1X B 


adjacent Mint midpoints »r* averaged to determine the values at aach jth 

station. Thus, 


♦ "3*3 


<3 • 2 . 3, . . N) (B5) 


(**b> 


"j-l*j-l*b,3-1 ♦ "j*j'b,j 


(3 • 2, 3, . . ., N) (B6) 


<**c> 


♦ "j*jrc,3 

3 * * 


(j » 2, 3, . . ., N) (B7) 


*3-3*3-^, 3-7 r c, 3-1 4 *3* j r b, j l c, j „ _ 

(ac b t c )^ - ~ (3 » 2, 3, ...» N) (B8) 


< x bb, j) 


*bb, j-J ♦ Ibb,j 


(3 • 2, 3 N) (B9) 


< I cc,3 l 


*cc,j-l * *cc, 3 


(j • 2, 3, . . N) (BIO) 


<*aa,j> 


*aa,3-1 ♦ J aa,3 


(j «• 2, 3, . . ., N) (BID 


. guation (B5) represents the lumped mass m-j at the jth station. For the end 
stations including the tip weight N 


(B)2a) 
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APPENDIX B 


*R*1 


«h*n w 
29 9 


(B’2b) 


*n*l*b,l 

^ • "IT" 


(Bl3t) 


(Wbijj+i 


"N«N*b,N ♦frb.N+l 

•♦■ 

29 9 


(Bl 3b) 


(•r c > \ 


*l*Uc,l 


C’Jl 4a) 




«N*N*C,N ^C^H+l 

4. 

29 9 


(Bl 4b) 


»»l*irb f irc,1 

(*br C )l 


(Bl5a) 


t**b C c) 


w N*N t b,N r c»N 

- 




^b,N+Uc,N+l 

9 




*bb ,1 



(B16«) 


*bb,N+l 


I bb,N*N 
+ 

2 


*n*n/*n 2 

I9 \i T 


2 \ 
r b,N/ 


Wfb,N+1 

4 

1 


(B’ 6b) 
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APPBtDIX B 


X CC,1 * 



4 ag \T r * 



(B17a) 


x ic,N%l *nW*N 2 , \ 

— ~ * 1? l u * '*•») * 1 — 


(B1 7b) 


x aa,l 




(B1 8a) 


txi 


x aa,H+l 


bb,N 


x cc^l)^l 


*H*N/ 2 

*irv tb ' H 


2 \ . #/ 2 
r C,N/ 4 ~\ r b,N+l 


r c,N+l) 


(B1 8b) 


Bqpiationc (B5) to (B18b) apply to the local axis system at each station. 
Because the equations of motion are written for the horisontal-vertical or x,Y- 
axea, the blade inertial properties must be rotated into the latter system using 
the pitch angle at each station defined from the twist 6 t and collective pitch 
angle e c . 


6j - 0 


(j » 1, 2, . . , , JHUB) (Bl 9a) 


® jHUB+1 * ®c 


<Bl 9b) 


® j+1 



(j - JHUB+1 , . . . , N) (B19c) 


where JHUB is an input quantity representing the number of nonfeathering hub 
segments. For general twist distributions, input values of twist are used, and 
the collective pitch angle is added to those values. The inertit representation 
in the X,Y,Z-axis system is given by the following relationships: 


(»r x ) ^ » (mr c )^ 

(mty)j • (mr c )j 


cos 6j + <mr b ).j 
sin 0j - (mr b )^ 


e<n 9j (j - 

cos 6j (j » 


1,2,.. ., Nil) 

1, 2, . . ., N+1) 


(B20) 
(B21 ) 
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*yy, j " *bb.j «* 2 ® j ♦ *cc,y • 1 » 2 ®j - 2<«b*c>j • 4 » e j °°« e j 

(i - 1, 2, . . .. M+l) (B22) 

l n ,j * *bb,s 811,2 ® j ♦ 1 cc, j cos 2 Oj ♦ 2(»r b r c ) ; j sin ® j cos ®j 

(j - 1, 2. . . ., M+l) (B23) 

*xf,i - < 4 bb, j - *cc,j> 8l « ®j cos ©j ♦ <«t>r c ) ^cos 2 ®j - sin 2 ®j) 

(j - 1, 2, . . H+l) (B24) 

Iu ( j * J-jj, j (j -1 # 2, . . .# H+l) (B25) 

* (Ibb,j - Icc,j) l 00 * 2 ®j ~ 8in2 ®j) ~ 4(*r b r c ).. sln ®j 008 e j 

(j » 1, 2 H+l) (B26) 

where the last ter* pertains to the so-called "tennis racquet* effect in the 
equation for torsional aoaent. 

The elastic properties are rotated into the X,Y,Z-axis system, using the 
total pitch angle "*t the midpoint of each segment obtained fro* averaging the 
values of pitch at the segnent endpoints. The shear center (elastic axis) off- 
s' ’£ fro* the reference axis (see fig. 2(b)) are found by 

i * i ®i 4 ®b r i 818 ®£ (i * !» 2* . • . , N) (B27) 

Sy # j * cos §j - S C/ £ sin 8j (i = 1, 2, . . N) (B28) 

The structural stiffness representation in the X,Y, Z-axis system is 
7 cos 2 8| sin 2 §i 

« + (i « 1, 2, . . ., H) (B29) 

<EI)yy,i (El) b, i <EI>c,i 
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APKMDXX B 


1 

<“>XX,i 


•in* 0j 
<*»b # i 


COB* 

<«>C,i 


(i • 1, 2, . . ., N) (B30) 


1 T 1 1 1 

— • - — sin 04 cos 04 

<«>sy,i [C*«b,i (K)c,iJ 


(i - 1, 2. . . ., N) (B31) 


The torsional stiffnesss (GJ) j is unchanged. 

At each station the centrifugal force and associated scwents created by 
the offset of the center of gravity from the reference axis are calculated by 
using the fonsulas 


11+1 

P j " 02 2 *k*k (j - n 2 N) (B32) 

k-j+1 

N+l 

F H»j * fl2 £ <>r X> k <j - 1. 2, . . ., H) (B33) 

k-j+1 


N+1 

p x,j ■ fl2 £ <**x> k *k (j - 1# 2, . . ») (B34) 

k»j+l 


N+1 

P y.j * fi2 S (wf y ) k* k <j I* 2, . . N) (B35) 

k-j+1 


where ft is the rotor rotational speed and z k is the radial distance to each 
station (fig. 1) 


*1 • 0.0 (D36a) 


i 

*j+i ■ £ * k 

k-1 


27 



APPENDIX C 


COMP UT ER PROGRAM LISTING 


The computer program listing is given in this appendix. 


*«3GRAN BL3 AML 1 iKPUTf OUTPUT f TAP E5»I**HTf TAP F6*0UTP0T# TAPE 7) BLOA 

C BLOA 

C THIS PROGRAM COMPUTES THE NATURAL FREOMFNCIFM AMO MOOC SHAPES BLOA 

C FOR A HELICOPTER POTO* BLADE. TEFTERINGf AMI ICULATEO* AMD BLOA 

C HINGELESS BLADES CAN BE TREATED. FURTHER INFORMATION BLOA 

C CONCERNING this PROGRAM NAT 8F FOOMO IN NASA TH 7B670 BLOA 

C BLOA 

CINNON /CO 1 * A/ jHHS # Nl*L0T»P3OTtITLEfl0l»NANE v 0AY»NPG BLOA 

I» JHOBI»RRPMaO)»RCOLLC3UZm>»INPUN»rfOOPLOT*N»TRSh»BLAOES BLOA 

COMMON/ CONC/ NTPLI 31I»EIB( 301 » UC( 301 .GJI 30) fTHD( 31 1 BLOA 

CQMMnN/CONB/lRCOLfkRCOL(3)f lRPN,CRP"(IO)f Z9AR(30)#EYEB(30)fSB<30)fBLOi 
i EYECOOI»SCI39lf SN # n)»RBUl>»RCOll>THI3I)»WTC30l BLOA 

C3M»ON/CaHO/CHATC5»*»'SDNNAmOO»3)»TPlNUOO»3), INODE ( 10Q# 3 > * BLOA 

1 «*3#HH5#CT(3n»5T(31)#l8»l<Tf IBS(3#10»3>f TBE(3f 10f3> BLOA 

COMMON/ k/ VMXnQ)fVOX(30>fV*Y(30)fF(30>f FTXC30)'FTY(30I#SX(30># BLOA 
I SY(30I# OMEGA?* FTnnOI BLOA 

COMMON/COMT / EYB(3n»EYCf31ltEYM3n»EMRB(31)»EMRCC3l)»EMRRm>» BLOA 
1 EMRX(31),EMRYtn>»ENBBIM3i>»ENPPM<31!»THHO<3i)»EMBPWt31> 8104 

COMMON /COME/ Z8(155lf Z«( 155 If ZQ< 155 If 7L( 155 >tZ$( 159 If ZY( 155 If BLOA 

i Z HUBS If ZD(155lf ZH(l55lf ZTU55I BLOA 

LOGICAL LOTtTRSN BLOA 

DIMENSION PP(2I# INOIf XOSOMI2lf$MZ(31lf SNZRYnU fSNZRY(31I BLOA 

I t SNR X ( ? 1 1 8LDA 

DATA CVRPS/O. 10*7198/ BLOA 

CALL PSEU03 BLOA 

CALL LEROY BLOA 

CALL D£TEC0A7| BLOA 

10 CALL 1NPT<plAST> BLOA 

MM 3*3 BLDA 

NH5*5 BLOA 

IF I TRSN I GO TO 20 BLOA 

NM3-2 BLOA 

NM5-* BLDA 

20 NOB«l BLOA 

.00 30 I* If 3 BLOA 

IN(I1«0 BLOA 

DO 30 'J‘1'10 BLC A 

00 30 X« If 3 BLOA 

1 SSI I f Jf K I *0 BLOA 

IBtUfJfKI »-l BLOA 

30 CONTINUE EcOA 

EN0X0»<PLASTACYRPSIPP2 BLDA 

00 25 i*lfI25f3I BLDA 


1 

2 

3 

* 

5 

6 
7 
B 
9 

10 

11 

12 
13 
l* 
19 
16 

17 

18 

19 

20 
21 
22 
23 
2 * 
29 
26 

27 

28 

29 

30 

31 

32 

33 
3 * 

35 

36 

37 

38 

39 
*0 
*1 
*2 
*3 


Z B 1 1 1 * 0.0 

zsm-o.o 

Z Y ( I I * 0.0 

zxm-o.o 

znm-o.o 

ZH(II*0.0 
Z0(II*0.0 
ZT ( 1 1*0.0 
ZKII- 0.0 
25 20(11*0.0 
Z B (1 1*1 «0' 

ZS( 321*1.0 
2 Y 1 125 1 *1 *0 
ZH 19* 1*1.0 
Z X (63 1 *1.0 

C COLLECTIVE ANGLE SWEEP * 
OQ 700 I$T a If IRCOL 


BLOA ** 
BLDA *5 
BLDA *6 
BLDA *7 
BLDA *8 
BLOA *9 
BLDA 50 
BLOA 51 
BLOA 52 
BLOA 53 
BLOA 5* 

ORlv tNAL PAGE IS l[li \\ 

OF POOR QUALITY \\ 

BIO* 59 

— — — t.Ot 60 

BID* 61 
BIO* 62 


28 



APPENDIX C 


€#*•**••#•*♦**♦**•***♦**•**•••♦***•**»•♦•*#*#•#********# 
C CALCULATE COEFFICIENTS DEPENDENT ON COLLECTIVE ANCLE » 
(**♦•*♦**•**#**•***»••***•*#*«******•**«»**♦*•#****•*♦•• 


90 


95 


ORIGINAL PAGE TS 
OF POOR QUALITY 


00 *0 I-KCB.Nl 

ZTX.THIII 

IMI.6T.JH<J«) ZTH*ZrH»XRCQL<I$TI 
XTII )•< INI ZTHl 

cmi-cosizrm 
SSTH-STII)P*Z 
CCTH.I.-SSTm 

sc?*«»sTm»cTin 

E*«li|.E«C(ll*CT(lN»«»UI*iT(ll 
EN 4 YIII-EPRCII >*STt II-EMRBI It*CTt I) 
E" 9 BWtt>«£YBtl)PCCTH»EVCII>*S'TH~Z.ftPE«RR(II*SCTH 
e it*Er)tn* 5 srH*Ercin*ccTM»z.o*EMi>in*scrH 
E 13 P 8 II l«IEVSt l l-EYCIll IPSCTH wEURRIIIPICCTh-STTHI 

TH‘o«n-ier 8 m-eYcu>»*<ccrH-$jTHi-*.o«F**t«m*scTs 

IFII.E 9 .Nll GO ro 80 
ZrH- 0 .*MTH(I>»THII»ll> 

IF ( I.GT • JHMBI ZTH>ZTH»«RCOlt ISTI 
STH- 5 INIZTHI 

cim-cosizthi 

SSTH«STH **2 

ccthm.-ssth 

SCTH«STH»CTH 

Sk(ii«scui»crH»$e<n*$rH 

SVtI)-S 9 tI>«CTH-SCII)«STH 

vnxu)>zb«r( n« (< i.o/ 6 ied i i-u. o/E tci I im* 5 cth 

VG<IIWURtI)*itSSTH/EIBtI>>»tCCTH/EICm>> 

V 6 r(ii*ZBk 4 m«((ccTH/Ei 8 (ni«<s$iHZFicmii 

CJNTIM'E 

SNRXINII- 0.0 

*HZIN 1 >« 0.0 

SNZRX(N 1 I» 0 . 

snz»Y<«i i«o. 

J-Nl t 
00 9 *> I» 1 »N 

J*J -1 

S' , «XI 4 >»S* 1 «XIKI*ENRXI»<I 

S«Z(JI«S"ZIX»*S 6 IX»*ZIX» 

5H2RXIJ>«SNZPXIk>*ENRJHK)PZI«I 

SNZRVt'l>«$*ZRYIK)«ENRYtK>*Zlft) 

•• 09 -JH 09 *! 


C ROTOR RPN $«E£P «- 

C*P** 


00 710 IB« 1 »IRPM 

C CALCULATE COEFFICIENTS DEPENDENT ON ROTOR R*N « 

c 

0 «E 6 AZ»CRPN(: 6 MPZ 

00 110 I» 1 *N 

FThI I) wONECAZ»S*RX(I> 

FI tl«ONEGAZ*SNZ(I ) 

FTXID-ONEGAZPSNZRXm 

FTrill-ONEGAZPSN/PUIl 

KAP$ 0 *I IEYEBII >»EYECII) )* 386 . 4 /WTPltI n«R 6 f t>*Rati>*RCtII*RCtt> 
WTIII- 0.0 

IFITRSNI wTtt) s Z 8 ARII) 7 tCJtI)*FtI)*KAP$Q) 

110 CONTINUE 

c •••••• ••*•••♦• ••*•••••»••••••••• 

C CALCULATE THE FIRST NATURAL FREOUENCY FOR EACH NODE TYPE • 

J«l 

C BYPASS COLLECTIVE ANO CYCLIC NOOES FOR NOPE THAN 2 BLADES 

IFIBLAOES.NE.Z.OI 4 »J 
00 ZOO I» J» 3 

IFIINIII.E 0 . 100 ) CO TO ZOO 


BLOA 

63 

BLOA 

66 

BLOA 

65 

BLOA 

66 

BLOA 

67 

BLOA 

66 

BLOA 

60 

BLOA 

70 

BLOA 

71 

BLOA 

72 

BLDA 

73 

BLOA 

7 * 

BLOA 

75 

BLOA 

76 

BLOA 

77 

BLDA 

76 

BLOA 

79 

BLOA 

BO 

BLOA 

BI 

BLOA 

87 

BLOA 

63 

BLOA 

66 

BLOA 

85 

BLOA 

86 

BLOA 

67 

BLOA 

66 

BLOA 

69 

BLOA 

90 

BLDA 

91 

BLOA 

92 

BLOA 

93 

BLOA 

96 

BLOA 

95 

BLOA 

96 

BLOA 

97 

BLDA 

98 

BLDA 

99 

BLOA 

100 

BLOA 

101 

BLOA 

102 

BLOA 

103 

BLDA 

106 

BLDA 

105 

BLDA 

106 

BLDA 

107 

—BLOA 

106 

BLDA 

109 

BLOA 

110 

BLOA 

111 

BLOA 

112 

BLOA 

113 

BLDA 

116 

BLDA 

115 

BLOA 

116 

BLDA 

117 

BLDA 

110 

BLOA 

119 

BLOA 

120 

BLDA 

121 

BLOA 

122 

BLDA 

123 

BLOA 

126 

BLDA 

125 

BLOA 

126 

BLOA 

127 

BLOA 

128 

BLOA 

129 

BLOA 

130 

BLOA 

131 
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*^n*U lM?0.U*CV*6S l •** 

*d$'l*t?»»00*0*CV*6$|99* 
cm CDFM U .f 4lSF.*2**0$0 **M* 

ic**«o 

WO 4*?10»*8$WM2>/*MiM 

IM04UQUT.1.0I GO TO 129 

I l O.U«Ml 

Mtn^Mn/ias^tn 
*p<2nppt2i/i9$ii'Wi 
go Yu wo 

1?9 MT»*»«M<?)*tv JS0«m~VQ$Q1<l t n«Q$PM2l 

If ^T«m< n«T4*P-kCSUlM?> WtQSOM?ll 
IF IT E $1 *1 T* 0 *C 0 1 1 GO n DO 
w? 

*9*P*M2nMT*4* 

W* C4U CPtHI»*F*tSC.»l»«Q$n*<?)»PP(2l» 

*P(?n*Pt:i/*951N<2l**tf *P 
GO T I :?* 

130 C0*T f NU*. 

huuiNtnu 

mu$?»ta»ii«i*in 

K-jcf*iMn 

SUNN* T I MHO* • W*F »T*U? 

*«Ef •S'INNaTINPODC * 11 / 2.0 
*QDF»N*P OF 

I F INODE .10.1001 CP TO 1*9 

C C 4LCUL 4t E H | CNF 4 H4TU94, ►HFQUfNClES FOO F4CM *00f TTPt ♦ 

li> 4QSu*ti n«3NNAt INODt* 11*1.01 
*3S( N|? l.SOHN*T|NODfc • 11*1.0? 

C4u ccmt»*F4i'f i^oso-ui.^un 
137 DENU*U**0S1N|1>**1£*P 
p£N0N«n*0$0MDI**lfIP 
DO 1*0 9 •NN33F • NPDt 

D* N"» L . q fc N3*U MAC 50*11* -SPUN *T 1**1 1*1.001 1 
1*0 DcMD*** OFHU^MFiJJKn-SOIMUn.n*! .ODD 
C 4 1 1 CPt F C 1 » ,t Al*c . »l»«QSOn:?» *P»<?)9 
f unc t •* * i n/oeNan 
rUNC*«*M?|/Ot*J»# 

? *!***• I <0>0"l >l**0.1*fl>NC«*t I |*OSC*i > >-3 4Fn**0.*>0-ClOS0»Ml 1 
1 tRF F »**U.SOI/lfl’NCt-FUHC«> ***•? 

If It «?«*?. *T.tNDXOI *0 TO HO 


0104 112 
0104 ill 
0104 134 
0104 111 
0104 114 
0104 117 
0104 DO 

old* ii9 
OLD" no 
0104 i .1 
01D4 142 
0104 141 
0L04 144 
Ot 04 141 
0104 146 
0104 14? 
OL 04 140 
OLD* 149 

old* no 

OL 04 111 
OL 04 112 
OL 04 111 
010* 11* 
*l 04 111 
OL04 116 
OL 04 11? 
old* ii* 
OLO* 119 
OL 0 4 160 
OLD* 161 
OLD* 16? 
OLD* 161 
OL 04 164 
0104 1 6 1 
OLD* 166 
OLO* 167 
0104 160 

sir* no 

OLD* 170 
OLD* 171 
OLD* 17? 
OLD* 171 
OLD* 174 
OLD* 171 


TFST«***Snt*t®4P-*QS34f?ll/!IQ$QOt2!l 
IMTFST.lT.O.OOll GO TO 160 

pptinppt?) 

»o>onih*)nso*i?i 

*TR4* 

GO TO 13? 

110 NODtnODf*! 

SONN*TI**lpE , l)*i||fUf 
X*EF»S*3NN*T|N0D£#U/2.O 
if noDtao.noi go io in 
GO TO 136 
1S9 CONTINUE 

N* ITF 16* 1H5 > 1 

18a FCRNATl//* l TE R *T ION $TPPP£1 f □« TYPE *#!?•* *OD£S 
IRLOUENCIEJ 6*3.*) 
no C 0**T l NUt 
IMt 3 • NPD* r 
10§II?T*10*n*"3Cl 
200 CONTINUE 

C CALCULATE and RRINT PUT NODE SHAPES ♦ 

CALL AHPlTD 
710 CONTINUE 


* L PA 176 
SlOA 177 
OLD* 170 
OLD* 179 
OLD* 100 
OLD* 101 
OLD* 102 
OLO* 1*3 
OLD* 104 
OLD* 101 
OLD* 186 
OLD* 187 
OLD* 100 
100 NATURAL FOLD* 109 
OLD* 190 
OLD* 191 
OLD* n? 
OLD* 193 
OLD* 19* 
OLO* 191 
OLD* 196 
OLD* 197 
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• .X* 

"too CONTINUE 

C MMF OUT S'JWUT OF NATURAL FREQUENCIES 

CALL SUNNY 

C HOT NATURAL FIEO. VS ROTOR IM • 

I M LOT I till RlWTIBlASTI 
60 TO 10 
ENO 


ORIGINAL PAGE IS 
1% POOR QUALITY 


•lor i?a 

6LCA 199 
6LDA 200 
6L0A 201 
BLOA 202 
6L0A 203 
OLD* 209 



SUBROUTINE anrltd 


ABPL 

1 

c< 



AWPl 

2 

c 

THIS SUBROUTINE CALCULATES AND RRINTS OUT MODE SNARES • 


AWPL 

3 

c< 



AWPL 

9 


CONNON /CONA/ JHU6*N1*L0T»R!3UT* ITlE (lH)#NA a E*0AY*NR6 


A9PL 

5 


1* JHUBl* RRRNI 1C 1 *RC(*LL ( 3) *Z ( 31 ) » INRUN* NODRlOT *N*TRSN»BL AOES 


AWXL 

6 


C0WH/COH5/l«COL#XRCOL(3l*lRPh*CRPW(l01*ZBAR(3OI#EYE6(3QI* 

SBC 30) 

*AWFL 

7 


1 ETEC(30>>SCI30)*SNI31)*AB(31l*RC(3ll*THI31)*VT(30) 


AWPL 

a 


C0NNQN/C0nD/CHAT(5*5),S0HNAT(100* 31* 1H«1 1W*,3I* INOOEt 100*31* 

AVIPL 

9 


1 nN3*NNS*CTI31)*STI31>*IB*IST,ISS(3»10*3).ISEI3* 10*31 


awpl 

10 


CORRON /CORE/ 29(1951# IXI155I# 20(1551* 21(1551# ZSU55I* 

ZY(155>AWPL 

11 


M1tl»l> Z0I155I. 2H( 1551* ZT 1 1591 


A«PL 

12 


C3NN0N /MINCES/ lch*lchri*lfh*lfhri ,chqff,fhoff*fcn,ffm 


AWPL 

13 


1 *R9HPUW#C , 3LPimLPH#LPHH*PH0FF*FPH* -OTS 


A9IXL 

19 


CONNON/CONT/ E YB( 31 ),EYC C 31 )»E YR ( 31 l»EHRB f 31) » ENRC ( 31 1* ENRR 1 31 ) » 

AWPl 

15 


1 EMt*(31),ENRV(31>»ENBBIH31)»EHRRN<31>*THHO(31t#ENBRMt31) 


AWPL 

16 


L06ICRL ROUT* TRSN* LOTS* NODRLOT * INRUN 


AWPL 

17 


DIMENSION VECISI* A(5)*I(3»31)»DNaT(5*5)» I RI V0TC5 I* IMM 101 


Awrt 

IS 


ORTR CVCRN/9.5A9297/ 


AWPL 

19 

c< 



AW7L 

20 

c 

RODE LOOP R*l FOX COLLECTIVE RODE ♦ 


AWPL 

21 

c 

*•2 FOR CYCLIC WOOES • 


AWPL 

22 

c 

*••3 FOR SCISSORS NODES • 


AWPL 

23 

c< 



AWPL 

29 


00 227 H* 1*3 


AWPL 

25 

c 

BTFRSS COLLECTIVE RNO CYCLIC NODES FOR A NUlTIBLADED ROTOR 


AWPL 

26 


IF(N.HE.3.AND.BIA0ES.NE.2.0> 60 TO 227 


AWPL 

27 


WOOEWO-O 


AWPL 

26 


IFlIBSt 1ST* IB, HI. EO. 01 60 TO 227 


AWPL 

29 

c< 



AWPL 

30 

c 

SHEER NATURAL FREQUENCIES STOREO IN SOHNAT • 


AWPl 

31 

c< 



AWPL 

32 


NRT - 1* IBE 1 1ST* IB*H) -IBS ( 1ST* IB*H| 


AWPl 

33 


00 223 NR*1*NRT 


AWPL 

39 


NRS • NR ♦IBS! 1ST* IB»H) -1 


AWPl 

35 


CALL COEF(N*. TRUE. *1»S0HN*T (NRS, N)*OUHYY» 


AWPL 

36 

c 



AWPl 

37 

c 

ITERATE TRIAL VECTOR FOR TIR DEFLECTIONS 


AWPl 

36 

c 



AWPL 

39 


00 9 I» 1 *HW5 


AWPl 

90 


DO 1 J" 1 * HW5 


AWPl 

91 


1 OHAT(I, JI-CNATI I* J) 


AWPl 

92 


A CONTINUE 


awpi 

93 


SQNNAT()>RS»Ni *CVCRN*S9RT (SONNAT (NRS»N)I 


XWP L 

99 


FNAT*SOHNAT(NRS,H) 


AWPL 

95 


CALL NATINV(5»NN5*CNAT*0»BBB»l,DETERN,irCAlE, IRIV0T.1HK) 


AWPl 

96 


VECdl-0.0 


AWPL 

97 


VEC(2)*0.0 


AWPl 

96 


DO 5 I ■ 3* HN5 


AWPl 

99 


5 VEC C I> *1 .0/l.E 50 


AWPL 

50 


ITftY-0 


AWPl 

51 


15 CONTINUE 


AWPl 

52 


ITRY*ITRY*l 


AWPl 

53 


ANAX*0.0 


awpi 

59 


00 25 I*1*NH5 


AWPl 

55 
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C 

c 

c 


Atmo.o 

00 20 4*WNH5 

20 m n*Atii*vecui*cwATu.4> 
iFtA0Sutm.0T.ASSi man 
25 CONTINUE 

DO 30 1»1»NN5 
iUMAtlt/llUV 
30 CONTINUE 

t F 1 1 TRY .EQ. A> GO TO 100 
00 AO 1*1 »HH5 
AO VECtl>*A( II/1.E50 
GO TO 15 


AHA L 
ANAL 5? 
ANAL 50 
ANAL 50 
ANAL 00 
AHAL 61 
AHAL 62 
AHAL 63 
AHAL 6A 
AHAL 65 
AHAL 66 
AHAL 6? 
AHAL 60 


CALCULATE HQOE SHAAES AND GENERALIZED INERTIA AHAL 

AHAL 

100 CONTINUE AHAL 

DO 60 I«l»H*5 AHAL 

VEC 1 1 1*0.0 AHAL 

00 50 4*1#HN5 AHAL 

VECClUVECtUOHATtmUAU) ANAL 

50 CONTINUE AHAL 

60 CONTINUE AHAL 

00 70 IM»NH5 AHAL 

At l>*At I)/AtNH5> AHAL 

70 CONTINUE AHAL 

l*Nl*l AHAL 

619*0*0 AHAL 

GI A»0.0 AHAL 

00 115 I*l»Nl AHAL 

1 *L-*1 AHAL 

G 180*281 nRAtll*ZBt3l*nRA<2UZBtl2A*I t*ZH«62*IIRAt3I AHAL 

G I AA* 2$tl)*A(lt*2$t3i + I}*At2! + 2Stl2* + I) + Z$t6? + Ik*At3! AHAL 

9 1 3#L > "0.0 AHAL 

Btl*l>*ZYtI>AAtmZYt31*I>*At2>*2Ytl29*mZYt62*l>RAt3> AHAL 

Bf 2»M*ZXtI>*Atl>*ZXt3l*I>*At2>*ZXtl?umiXt62*l)RAt3) AHAL 

IM.NOT.TRSN> GO TO 113 ANAL 

G16B*GlBB + ZBf93M)*AUI ANAL 

GIAR»GlAR«ZSC93*n*AU> AHAL 

B t l»L)*8tl»L>+ZYt93+I)+At9| AHAL 

9t2'M*6f2»L>*ZXt93*l)*AU> AHAL 

Bt3tL>"ZHtI>RAtlUZMt31*l>*At?>*ZHtl2A*I>*ZM<62*Il*At3>*ZHt93*I>R AHAL 

1 At A I AHAL 

113 CONTINUE AHAL 

G I B-GI B + EHBBNt L I + GIBB**? ANAL 


69 

70 

71 

72 

73 
7* 

75 

76 

77 

78 

79 
BO 
81 
62 
63 
BA 

65 
St> 
87 

66 

89 

90 

91 

92 

93 
99 

95 
9t 
97 

96 
99 


GI A»GIA^ENAANtL>RGIAAR*2 
115 CONTINUE 
AB SCL*0 . 

00 120 I*1»HH3 
00 120 4 • 1 • HI 
ABS6* ABS IBtlfJI) 
IEtABSCL.GT.ABSB> GO TO 120 
SCALE-Btl# J) 

A BSC 1 • ABSB 
120 CONTINUE 

SCALE-1. /SCALE 
GI8»G16*SCALE**2 

gia«gia*scaler*z 

DO 150 1*1 » HH3 
DO 150 J-lfNl 
150 Bt I» 4 1 -Bt I.J)*SCALE 
GlX-0.0 
GIY-O.O 
GIT-0.0 
DO 317 J-UNl 
GIY-GlY^SHt J)«B(hJ)M2 
GiX-GIXt$NU)RB(2» J>R92 


AHAL 100 
ANAL 101 
AHAL 102 
ANAL 103 
AHAL 109 
AHAL 105 
AHAL 106 
AHAL 107 
AHAL 106 
ANAL 109 
AHAL 110 
AHAL 111 
AHAL 112 
AHAL 113 
AHAL 1 1 A 
AHAL 115 
AHAL 116 
AHAL 117 
AHAL 116 
AHAL 119 
AHAL 120 
AHAL 121 



APPENDIX C 


C 


c 


c 


ORIGINAL PAGE £ 
OF POOR QUAIOT 


GIH.tT.6lT) m*t*PS>*)»? 
GIT.GT.GIH) IPLNlHtS »*)•* 


GIT-GIT^ETRU)* 8(3*J )**2 
317 CWtHUE 

GIV*G18*GIY 
GI^GIP+GIX 
IPLMNPS*M-l 
I F (GIH.GT.G l V . AKO* 

IF (GIT.GT.GI V . ANO. 

FREQPR»FNAT 

IF(Rft?H( tel.NE.O.O) FREOPR*FNAT/RRPH( 18 ) 
emTE NUHBEA CF HOOFS 

IN10E(NPS»8I-0 

Kft*JHOBl 

IF(JHUB.FO.O) **-2 
DO 160 l*PX*N 
PR00«3(R, n*B(K«I*l> 

I F (PROO.GE .0*0) GO TO 160 
INODE (NP$*8) • INQDE(NPS*N ) +1 
160 CONTINUE 

IF (LOTS) GO TO 190 
IF (.NOT. POUT) GO TO 225 

IF(RCOLL(IST).NE.COLPUN.OR.RRPN(IB).NE.RP*PUNI GO TO 225 
PLOT HOOF SHAPE IF REOUESTEO 

IFIHOOPLOT) CALL POOSHAPf 8*RC0LL ( IST)*RRPH(I8 I * FR£OPR*H*NP# HPT) 
190 CONTINUE 

PRINT OUT NODE SHAPE 
NPS*HPGM 

WR ITE (6*901 ) NPG*DAY#NAHE#(ITL£(4)#J»1*10I 


AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 

AHPL 


901 F0«HAmHl*l8X*4HPAGE*l3*25X*34HHYRLESTA0 ANALYSIS NASA TH 7B670 AHPL 

1 *25X»AlQ»//10X*5HCASE * 16. 23X, 10A6/ ) AHPL 

HR ITE (6*904) RCOLL ( 1ST) *RRPH( IB) AHPL 

904 F0RHAT(37X*F6.2»2X*23HDEGREE$ ROOT COLLECTIVE* 10X*F7.1#2X» 16HR0T0RAHPL 
1 SPEEO* RPR) AHPL 

I F ( ».FQ*l) MR ITE ( 6* 902 ) FNAT AHPL 

902 FORHAT ( 36X* 27HC0LLECTI VE HOOE OF BLADE AT »F9.2*4H CPH ) AHPL 

I F ( N.E0.2I WRITE (6*903) FNAT AHPL 

903 FORHAT ( 38X* 23HCYC*. IC HOOE OF BLADE AT *F9.?*4H CPH ) AHPL 

lFfH.EQ.3) WR ITE (6*916) FNAT AHPL 

916 FORHAT ( 38X* 25HSC ISSORS HOOE OF BLADE AT * F9.2*4H CPH ) AHPL 

IF(RRPM( IB). HE. 0*0) WRITE(6* 920) FREQPR AHPL 

920 F0RHAT(IH**76X*3H0R *F?.3*8H PER REV) AHPL 

WR ITE (6*909) AHPL 

909 FORHAT ( /10X* 9H8L ADE STA*BX* 11H0E FLECTIONS * 19X * 7MHGHEF V S *20X* AHPL 

I 12HSHEAR FORCES *12X*5HTWIST *6X*6MT0RQUE /13X*2HIN, 16X*2HIH * AHPL 


122 

123 

124 

125 

126 
127 
12 B 

129 

130 

131 

132 

133 

134 

135 

136 

137 
130 

139 

140 

141 

142 

143 

144 

145 

146 

147 
140 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 


2 2 5X# 5H IN-LB* 26X* 5HLB * 15X* 3HRA0* 7X* 5HIN-LB* /26X* AHPL 

3 4HVERT* 4X* 5HH0R IZ*13X*4HBEAH*8X* SHChORO* 13X*4H8E AH*6X* 5HCH0R0 / AHPL 

4 2X*16(BH*P44RAR* )) AHPL 

L*Nl4l AHPL 

00 220 J" 1 * Nl AHPL 

1*1*1 AHPL 

BH«(2HIL)PA(1)4ZH(31^L)9A(2)*2HU244L)47H(62*L)9AI3))9SCALE AHPL 

0-(20(L)4A(l)4Z9(3lPL)PA(2)^ZO(124*L)420(62tL)4A(3))*SCALE AHPL 

EL-(ZL(L)AAa)4ZLl3l4L)4A(2)4ZL(1244L)4ZL(624L)4A(3))4SCALE AHPL 

0E-(Z0(L)»A(l)*Z0(3L4L)*A(2)4ZDtl244L)4ZD(62H)*A(3))*SCALE AHPL 

IF(.NOT.TRSN) GO TO 235 AHPL 

T ■<ZT(L)PA(1)4ZT(314L)*AC2)4ZTU244L)4ZT<62U)9A(3)4ZT< 93H ) *A( 4AHPL 

1 ) )*SC ALE AHPL 

BH»BN4ZH(93H)*A(4)*SCALE AHPL 

0 «0 +Z0(93*L)*A(4)*SCALE AHPL 

EL»EL4ZL(93+L)*A(4)*SCALE AHPL 

0E*0E4Z0(934L)9A(4) RSCALE AHPL 

235 CONTINUE AHPL 

IF (PHOFF .HE .0*0 .AND. J.LE.LPH ) T*0.0 AHPL 

QBN>SH*CT(J)+Q*STU) AHPL 

OCO-OACT(J)*BH*ST(J) AHPL 

FBH-ELACT(4)40EaST(J) AHPL 


166 
167 
16 B 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 
160 
181 
182 

183 

184 

185 

186 
167 


33 
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APPENDIX C 


FC0»DE*CT(4)-EL*$T(4) 

u*m<6»<un 4* iui# su.ji* M2»4i* os**» ocd* mu* 
U1 FQRNAY(6V*t3»F8+2*FlA*3»F6*3tFl9,0»F12»0»Fl6«0»FLl*0) 
1FI TtSN I WR1T€<6*912) 6(3*41* T 
112 FORMAT tlN* f 971 *F 13. 1*F13.0) 

220 CONTINUE 

WR1TH6»91A) 

91A F0RNAm6(8NA***A**A>) 

WRtTF(6.A0?l 6ltf»6lH*GlT*(VECtn*l«l*MN5) 

A07 Fa**AT«/10X*35HGENERALT2EO INERTIA UN-LNF-SECA*2 )* 
l 5*AA»5HH0RH*F9.5 # 5X,5HTWlST*FB.3*7*10X*2AHMn)NDAaY 
?C«*3X*5E12.5> 

223 COhTJNUF 

IFUNUT.IMUNI C>0 TO 223 
•HNCH OUT «nOF SHARE 

lF(Nf)0EN0.LT.6) CALI CAROS (N*N0DEN0*MRS*6* 1 > 

{23 CONTINUE 
{27 CONTINUE 
RETURN 
€N0 


ANRL 166 
ECO AML li« 

AML 190 
ANAL 191 
ANAL 192 
ANAL 193 
ANAL 19A 
ANAL 193 
ANAL 196 
3X*AHVERT*F6.ANRL 197 
COMO IT ION CHCANAL 199 
ANAL 199 
ANAL 200 
ANAL 201 
ANAL 202 
ANAL 203 
ANAL 20* 
ANAL 203 
ANAL 206 
ANAL 207 


subroutine C ARO$(N*HOOENO#NB$* 6* I F 18ST) CARO 1 

CARO 2 

This SUBROUTINE AUNCHES OUT NODE SHARES • CARO 3 

CARO A 

COMON /coma/ 4Hi'B»Ml*LOT»ROOT»lTLEUO>*NAME»OAY»NRa CARO 3 

l*4HUBl*RRRNI 10) *RCOLL(3l*2tU)* INRUN, N00RLOT*N*TRSN*6LA0ES CARO 6 

COHNPN/ CONC / WTRL 1 31) *E I6( 30) »E t CI30) »<>4 1 30 ) »ThO( 31 ) CARD 7 

CONAON/CON9/1RCOL* RRCOL <31* 1RAN»CRAN| 10) » ZB At 130 1* I YEB 1 30 1 * SBI 30 1 *C ARO 8 

1 E TEC (30)*$C(30)*$Nt3i)*R8(31)*RCt311*THt31)*ltT(30) CARD 9 

C ONNOW /CONO/C NAT ( 3* 3) fSONMAT ( 100* 3) • t ALMf 100* 31* INOOE ( 100* 3 > * CARO 10 

1 NR3 # RH5»CT<31)»S1 (31) *18* 1ST* I BS<3* 10*31 ilfttt 3*10*3) CARO 11 

CONNON /HtNGES/ L CH» LCHA 1* LFH* LFHA1 * CHOF F * FHOFF * FCH* FFH CARO 12 

1 *RANAUN*COLAUN*lAH*LAHAl»AHOFF*FAM* LOTS CARO 13 

0 1 NENS ION 4(2 l)*AHQUX(21)*RHOtYT(211»B( 3* 31 )*AUM(22* 31*0(6*3*6) CARO 1A 

OATA W* RHOITX* RHH lYY/63*0 .0/ CARO 13 

IFUFIRST.GT.O) CO TO 33 CARO 16 

CARO 17 

RECAST BLADE RUNNtNG WEIGHT ANO INERTIA DATA IN10 20 EQUAL CARO 16 

SEGMENTS CARO 19 


N8I1**7*D NA^fetl 1TLE( I ) * t • L * 1 0 ) 
l F 1 '*».>iU6*AY»10AM 

W<n)»wT*L(AU 

0* 11**121 )*«fr>L(M)»AB(Nl)*«2 
5HllYT(21)»4T» > l(M)AaC(NlMA2 

3EL9«2(M)/20«0 
J» 1 

fc(l)*0.0 
RhOIXK 1) *0*0 
jiyy \ 

IS I«l*20 
* TATP«TRDFI» 

5 CONTINUE 

1 F (7CJOI .J) • S T AT r 1 GU TO 10 
W( l)»wm ♦JTRL < J) •, 8AM J ) 

<rt HX*U )»9HCll«M II + (tYEM4>*NTRL(4l*R*(4l**?/366.6)*2BAft(4) 
RH JlYYI l)«JHCIYTI 1 ) ♦€ EYEC (J3 ♦UT°L(J 1*PC IJ)**?/T86*A)AZ&ARiJ ) 
J • J ♦ 1 

I M M J I O Q.SfATJ) GO TO 13 

g? n 3 

10 h e AC - ^ T AT L-7 ( U ) 

W(I)*ifUI*rfiAC*4TAL(4) 

R HO I Yin I»9HU1U< n*MAC*<fYFMj)*wT*L( J) f RR(J)AA2/3B6#A ) 


CARO 20 
CARO 21 
CARO 22 
CARO 23 
CARO 2 A 
CARO 25 
CARO 26 
CARD 27 
CARO 2 8 
CAPO 29 
CARO 30 
CARO 31 
CARO 32 
CARO 33 
CARO 3 A 
CARO 35 
CARO 36 
CARO 37 
CARD 38 
CARO 39 
CARD AO 
CARO A 1 
CARD A 2 


34 
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APPENDIX C 


ftHaiTV(l)-ftH3IY7(Il*Flue*(EVEC(4)*WTPl(J)*ftC(J)**2/3e6.0 
w(i.i)*fiuc*wtpl< ji 

ft HQ t XX ( I*l)*FftAC* (ETE8( 4)*bTPl (4 ) *1(8(4) ft*2/3ft6»4 1 
ftHJIYY( 1*1 )*FftAC*(EYEC(4)+WTPL (4 )*ftC (4)**2/386«* t 
4*4*1 


19 CONTINUE 

OQ 20 1*1.20 

W(t)*W(I)/OElft 

RH01XXU»*»HOtXXm/0Etft 

20 «HQIYY(I)*ftHOtYY(n/OEl« 

WRITE (7.301 (Wilt. 1*1. 21) 

WRITE (7.30) (ft HO IXX(t). 1*1.211 
MR ITE (7.30) (RH0IYY(I).I*1.21) 

30 FORNftT) 7F10.5) 

RETURN 
35 CONTINUE 


./v.iNAL PAGE IB 
POOR QUALITY 


RECAST NOOE SHARE INTO 20 EQUAL SEGMENTS 
NOOENO*NODENO*1 

IF (RCOLK 1ST ) •NE.COlftUN.OR.RRPM IB) .NE .RRNRUN) 60 TO 702 

00 AO «*1#3 

RUN(l.K|*B(K.l) 

AO RUN(21.K)*B(K.N1) 

0ELft*Z(Nl)/20.0 

4*1 

00 95 1*2.20 
STATO*( I-1)*0ELR 
A5 CONTINUE 

1F(Z(4*1).GT.ST*T0> GO TO 50 
4*4*1 

1F(Z(4).LT.STAT3) GO TO A5 
00 A7 K*1.3 
A7 RUN(I.K)*B(K.J) 

GO TO 55 

50 FRAC*(STAT0-Z(J))/ZBAR(4) 

00 52 K-1,3 

52 PUN(t.K)>B(K,J).FRAC*(B(K.4«l)-B(K. J)) 

55 CONTINUE 
00 60 1*1.21 
PUN(Itl)-PUN(l.l)/12.0 
PUN(1.2)*PUN(I. 21/12.0 
60 PUN(1»3)*PUN(I.3)«57.295B 


fnat*SgnnaT(n»s#n> 

IF(RRPH( IB) .NE.O.O) FNAT*FNaT/RRRN(IR) 

PUN( 22* 1 ) *FN*T 
I E (ft. E4 • 1MUN(22,2»* 1. 

IF(N.EO. 2 1 PUM 2? .2 ) «-l . 

1F(«.E0. 3 ) PUN ( 22 .2 )*0. 

*UN( 22. 3 ) * . 02 

IF«ABS(Pl'N(22.l)-1.0).LE. 0.05. AND. «.E0.») GO TO 200 
GO TO AGO 

200 00 A50 AN. 1. 22 
PUNtftN. 3>*0.0 
A50 CONTINUE 
A60 CONTINUE 

ISW*PUN(22,2) 

00 600 KKft*1.21.2 

KKKPl.ftKK .1 

NR ITE (7. 27 » (PUN(RKN.I).I«1.3) » (PUN (KARRI* II. 1*1 . 31 .NAME. 
1 NODENO.ISN 
600 CONTINUE 
RETURN 

702 CONTINUE 

C PUNCH CYCLIC DETUNING CARO FOR EACH NOOE SHAPE 


CARD 

A3 

CARD 

A* 

CARO 

A3 

CARO 

A6 

CARO 

A7 

CARO 

AS 

CARO 

A9 

CARO 

30 

CARO 

51 

CARD 

52 

CARD 

53 

CARO 

5A 

CARD 

55 

CARO 

56 

CARO 

57 

CARO 

5S 

CARD 

59 

CARO 

60 

CARD 

61 

CARO 

62 

CARO 

63 

CARO 

6A 

CARO 

65 

CARO 

66 

CARO 

67 

CARO 

66 

CARO 

69 

CARO 

70 

CARO 

71 

CARO 

72 

CARD 

73 

CARO 

7A 

CARO 

75 

CARO 

76 

CARO 

77 

CARO 

78 

CARD 

79 

CARO 

SO 

CARO 

81 

CARO 

82 

CARD 

83 

CARD 

8 A 

CARO 

85 

CARD 

86 


CARO 

37 

CARO 

8 8 

CARO 

89 

CARD 

90 

CARD 

91 

CARO 

9? 

CARO 

93 

CARD 

9* 

CAPO 

95 

CARO 

96 

CARO 

97 

CARO 

98 

CARD 

99 

CARD 

100 

CARD 

101 

CARD 

102 

CARO 

103 

CARD 

10 A 

CARD 

105 

CARD 

106 

CARO 

107 

CARD 

108 
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APPENDIX c 


IF(ia.EQ.l.AND.I$T.EQ.l)0(l.H,NODENn)>SQHNATfNRS.H> 
IMI9.E0. 1. AND. IST.EO. 3)D( 2.H. HODEND-SINNATINRS.H) 
IHI8.EO.3.AND.IST.E0.1)0(3»N,HODENO>-*ONNAT(NPS.H) 
IF(IB«EQ«3.AN0. tST<EO»3)D(A.H.HOOENP)- r OHNAT(NP$.H) 

27 F0RHATI6F10.6.1 OX. 16.212) 

7 CONTINUE 
L- 1 

IF(BLA0£S«GT.2.0) 1-3 
00 8 IK-1.3 
00 8 I J-1.6 

01 5* IN, I J ) -(RCOLL ( 3 >-RCOLL ( 1 ) I *. 5 
0(6.IN,IJ)-(RRPN(3)— RRRH(1))*.5 
IFIIN.EO.ll ISW-1 
IFIIN.E0.2i !SW — 1 
IFdN.E0.3I ISW >0 

IF (IB.EQ.3.AND. IS T. EO . 3 . AND .NODENO.EO .6. ANO.N . EO . 3 ) 

1UR ITE ( 7.28 ) (0(l.IH.IJ).I-1.6)»RC0LL(2)«RRPH(2).NAHE.IJ.ISW 

28 F0RNATI6F 10.2. 2F3.0. 16.212) 

8 CONTINUE 
RETURN 
ENO 


CARO 109 
CARD 110 
CARO 111 
CARO 112 
CARO 113 
CARO 1 1 A 
CARD 113 
CARO 116 
CARO 117 
CARD 118 
CARO 119 
CARO 120 
CARO 121 
CARO 122 
CARO 123 
CARO 12 A 
CARO 12* 
CARD 126 
CARD 127 
CARD 128 
CARO 129 


SUBROUTINE COEFI IA , 


OET* INAk. SOSON, Rp) 


THIS SUBROUTINE CALCULATES THE DEFLECTION OF EACH STATION 
AS A FUNCTION OF DEFLECTIONS AT THE ROTOR BLADE TIF FOR 

IA-1 FOR COLLECTIVE NODES 
IA-2 FOR CYCLIC NODES 
I A-3 FOR SCISSORS NODES 

INAX-NO OF FREQUENCIES TO SC CALCULATED 
SQS0HI1 TO INAX) CONTAINES SQUARES OF FREQUENCIES 

DET-.TRUE. USED TO FIND NODE SHARE FOR KNOWN NATURAL FREQ. 

OET-. FALSE. USED TO FIND THE OETERNINANTS OF THE BOUNDARY 
CONDITION NATRICIES 


COHHON /CONA/ JHUB.N1.L0T. POUT .ITLEI10i.NANE.0AY.NRG 


1. JHU81,RRPN(10).RCOLL(3).2(31).INRUN,NOORLOT.N.TRSF .BLADES 
C0NH0N/C0N9 /IRCOL. XRCOL ( 3 ) . IRRN. CRRN ( 10) . 2BAR ( 30 ) . E YEB< 30). SBC 301 
1 E YEC ( 30) .SC ( 30) . SN ( 31 | .RB ( 31 ) .RC (31 i.THI 31 ). WT( 30) 

CONHON /COHO /CHAT (3. 3).S0NNAT( 100. 3).IRLN(100.3). INODE (100. 3 ) . 

1 NH3.NH5.CT( 311. ST(3\)» IB. 1ST. IBSO. 10.3). IBE(3. 10.3) 

CONHON/ H / VNX(30).V0X( 30). VNY (30).F(30).FTX(30).FTV(30).SX(30). 

1 SY(30).0HEGA2,FTH(30> 

CONHON/ COHJ/ KIP.NIR.KOR.HOR.KBETA.KRSI.KT.KC 
CONHON /HINGES/ LCH. LCHR1, LFH. LFHR1 . CHQFF.FHOFF.FCH.FFH 
1 .RRHRUN,C0LRUN,LRH,LRHR1»RH0FF.FRH, LOTS 

CUNHON/COHT / EYB(31).EYC(31).EYR( 31 ) . EHRB( 31 ) .ENRC ( 31 ). EHRR ( 31 ) . 

1 ENRX ( 31 ) .EHRY ( 31 ), EHBBWI 31 ) . EHRRW( 31 ).THHO( 31 ). EHBRWOl ) 

CONHON /CONE/ 28(155). ZX( 155) . 20(155). 2LI155). ZSU.5). ZY(155 
*» ZN( 155 I. 20(155). 2H(155). 2T(155) 

REAL KIP.HIR.KOR.HOP.KBETA.KPSI.KT.KC 
LOGICAL OET. TRSN 

D INENS ION SQS0H(1)»PP(2>»IPIV0T(5).IWK(10) 

C CALCULATE DEFLECTION COEFFICIENTS * 

00 220 II-l.IHAX 
SOHS-SQSOH(II) 

SONSY- SONS PON EG A2 
ZN(1)-$0HSY*EHBBW(N1I 
ZQ(1)-S0HS*EHBPH(N1) 


COEF 

1 

COEF 

2 

COEF 

3 

COEF 

4 

COEF 

5 

COEF 

6 

COEF 

7 

COEF 

6 

COEF 

9 

COEF 

10 

COEF 

u 

COEF 

12 

COEF 

13 

COEF 

14 

COEF 

15 

COEF 

16 

COEF 

17 

COEF 

IS 

COEF 

19 

i#COEF 

20 

COEF 

21 

COEF 

22 

COFF 

23 

COEF 

24 

COEF 

25 

COEF 

26 

COEF 

27 

COEF 

28 

COEF 

29 

COEF 

30 

>)COEF 

31 

COEF 

32 

COEF 

33 

COEF 

34 

COEF 

35 

COEF 

3b 

COEF 

37 

COEF 

36 

COEF 

39 

COEF 

40 

COEF 

41 
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APPENDIX C 


ORIGINAL PAGE IS 
OF POOR QUALITY 


ZH<32)*$0NSY*EN8PV(N1I 
2a(32»on»$*EN»pw<tin 
ZDl63>»SUN$Y*$N<Nl) 

2T<63)*SUHSY*ENRV(M) 

ZM94I— UPEGA2*2(m*£NR* .Nil 
20(94)— nNeGA2*Z(Nn*CHRY(Nl> 

ZLt 94) •$0*$*F**V (Nl) 

ZD I 94|.$0iSY*E*RY(Nll 
zn 94)«SON$*EYR<Nl)*UNF6A2*THM0(Nll 
2L(125>*$H(Nl)*SOHS 
ZT(125)*SON$*EMftX(Nl) 

00 135 j* 1 • 5 

tF(.NQT.TRSN .AND. J.E0.4) GO TO 135 
*«N1 

L1OO1-30 
00 X 30 I*2*Nl 
1 2 *1 X 
L 1 * L 1 ♦ 1 
N«*-X 

ZQVTU*ZBA«(N)/2.0 
2Sa$X«Z8AR(H>**2/6.0 
V*HEAR-2L(L2)-F(*)*ZB(L2) 

H$HEAR«ZD(L2)-F(H)*Z$(L2I 

OH*WT(N)MH$HEAr*$Y(«| - VSHEAR*SX(N| ♦ FTY CHI *2$ t 12 1 ♦ FTX(H)* 

X 28(12) ♦ ZH(L2I*$X(N>9FTH(N> - 2TU2II 
08— ZOVTO*(VNY<N)*VSHEAR*VHX(H>*HSHEAR) - VNY(N)*2H(L2l - VHX(HI* COEF 
X 20*121 - ZM(L2>*(f rY(*)«VttX(*)*FTX(N)*VNYf«)-2OVT0*FTH(*)*VHY(NnCOEF 
OS— 20VTUMVNX(N)*VSHEAR*VOX(N)*H$HFAR| - V*X(N)*2ML2) - VOXXIU* COEF 
X 20(12) - 2H(l2)*tf TY(N)tVQX(N)*FTXt»)*VHX(H)~ZOVrO*F?H(ft)*VftXCN))COEF 
DY*Z$QSX*(VHY(8)*V$HEAR*VHXCH)*HSH6ARJ - 28 AR (R ) *28 ( L2 I ♦ VHY(N)* COEF 

1 2 OVTO* ZN ( l ? I ♦ VhX(HI*2QVT0*20(L2l ♦ 2Hf 12 )* (FT Y (IO*ZOVTO*VNXOO COEF 

2 ♦FTX(N)*73VT0*VPY(H)-FTH(H)*Z$QSX*VNY(N>) ♦ SX(«|*OH COEF 

DX*ZS0SX*(V9X(N)*H?H£AR*VHX(N)*V$H£AR) - f 8 AR (H| *Z$(L2 I ♦ V«X(N>* COEF 

1 ZflVTO* ZM t L 2 ) ♦ VQX(*I*ZOVTO*ZQ(L2) ♦ ZM( l? I ♦ (FTY( N >*ZOVTO*VOX ( *> COEF 


COEF 

COEF 

COEF 

COEF 

COEF 

COEF 

COEF 

COEF 

COEF 

COEF 

COEF 

COEF 

COEF 

COEF 

COEF 

COEF 

COEF 

COEF 

COEF 

COEF 

COEF 

COEF 

COEF 

COEF 

COEF 


2 ♦ FTX('*)*ZOVTO*VPX(N)-FTH(N)*Z$Q$X*VNX(N> ) - .>Y(H)*DH 
tHBl >t e .E0.2.0l GO TO X20 

IFfFHw, F.EO.O.O ♦ AND .CHOFF . £0. 0. 0 ) GO TO 120 
iFtZOW.LT.FHQFF) DB — 2B(L2) 

1MZCR) .LT.FHOFF) DY — ZY(L2) 

IF(7(8).LT. CHOFF) D$ — ZSU2) 

IF(Z(N).LT. CHOFF) OX— Z X ( L2 1 
120 Z3(L1)»28(12)*0B 
Z$ (LI )*ZS (L2I+0S 
2Y(L1)«ZY(L2)*0Y 


COEF 

COIF 

COEF 

COEF 

COEF 

COEF 

COEF 

COEF 

COEF 

COEF 


*2 

43 

44 

45 

46 

47 

48 

49 

50 
5 X 

52 

53 

54 

55 

56 

57 
5 R 

59 

60 
61 
62 

63 

64 

65 

66 
67 
66 

69 

70 

71 

72 
7j 

74 

75 

76 

77 
76 
79 
60 
81 
62 
83 

64 

65 


Z X (L l ) > 2X 1 1 2 ) *D X 
2H(L1)-2H(L2)*0H 

U (Ll )*Zt(l2>*$0N$* (ERR X (HI *ZH< LX )*$N(N>*2Y( LI H 
ZD(UI*Z0(l2>*S0H$Y*(£*RY<N>*ZM<ll>*?AC«O*Z)r<ll> I 

ZMLll*ZN(L2l*F(l1>*DY*ZBAR(N)*2L(l2l*$ONSY*(EN8PW(H>*ZS(ll>*EHBBN(COEr 

1 ft)*2B(ll) )-DNEGA2*Z(R) #ENRX(N)*2H(ll) COEF 

Z0(ll)*ZQ<l2)*F(M)*0X*Z8AR(H)*ZD(l2)*S0N$*(FNPPW(N>*ZS<ll)*E«BP*(HC0EF 

2) *26(11 ) >-3HEGA2*ZM)*EflRYt*t*2H(ll) COEF 

ZT(lll*ZT(L2l*FTM(NI *OY*SON$Y*ENRY( Nl *ZX* ll )+SQNS*ENRX (NI+ZYILL) COEF 

2 ♦(SOHS*FTR (*U*0«EGA2*THH0(R|)*ZHCtll 
130 CONTINUE 
135 CONTINUE 

: CALCULATE ROOT CONDITIONS PER INDIVIDUAL TIP DFFLFCTI0N5 
L-Nl-31 

DO 200 J-1#NN5 
1*101 

IF(.NOT.TRSN) GO TO 193 
: TORSION CONPONENT SOUND AH CONDITION 

CNAT<5# J)«ZT(L-lPH)-FPH*(2T(L-LPH)-7T(L-LPHPll)-KC*(2H(L-LPH)-FPM*COEF 104 
lC2H(l-LPH)-ZK(L-lPHPl)n COEF 105 

GO TO 197 COEF 106 

193 Z F (J • EO* 4 ) 1*101 COEF 107 


COEF 86 
COE r 67 
COEF 66 
COEF 69 

90 

91 

92 

93 

94 

COEF 95 
COEF 96 
COEF 97 
COEF 98 
COEF 99 
COEF 100 
COEF 101 
COEF 102 
COEF 103 
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197 CONTINUE 

IMU-2) 6 0.5D0#600 
600 CONTINUE 

c boundary conditions for collective nodes 

CPAT<2#4) • 2MU 
C^AT ( 3* J ) • ZSill 
CNATt6#J|.ZQ<l)-KT*ZS(U 
GO TO 200 
500 CONTINUE 

C 80UN0ARY CONDITIONS FOR CYCLIC NODES 

• mu 

CNATC2#4I«20(U-|KiP-S0NS#NlP>A2XU) 

CPAT<3»JWPUI-KBETAtZB(l) 

CPAT<6.3) • ZMtl 
GO TO 200 
600 CONTINUE 

C BOUNDARY CONDITIONS FOR SCISSORS NODE 

C«Am#J)*ir,:-lFH>-FFHAIZYtl-l FH>-ZY(l-LFHPin 
CNATI2#J>«2*U-ICH)-FCHMZXU-LCH)-M<t-LCHM H 
CNATOtJ) • Z9U) 

CNATIA.JJ • ZMl) 

IFITRSN * AND * Bl ATE5 *£0* 2*0) 

1 CBATt5#J>«ZHU-lPH>-F»H*<IHU-lPHI-ZH(l~LPNPLn 

I F ( FH3FF *07*0*0) CNATn#J>«ZN(L-LFH>-FFM*(ZNU-LFH|-ZHlL-LFHPin 
1 -K8ETAAlZB(l-LFH|-FFHAUB<L-lFHl-2ML-tFHPnM 
I F (CMOFF . GT ♦ 0* 0 l CPAT<6#J)«ZOtl-lCM|-FeMMZO<L-tCH|-ZQU-lCMPl > ) 
l •<PSXRIZSCL-LCPl-FCHA(ZSU-lCHI-ZS(l-tCH»l)M 
200 CONTINUE 

I F COE T I RETURN 

CAU HATINV(5*NN6»CNAT»0>BBft»0»DETFRN»lSCAlE»IPl VOT # INK ) 

P* II II «0E TER** 1Q.Q*M 100*1 SCALE) 

220 CONTINUE 
RETURN 
END 


COE F 106 
COIF 109 

coif no 

COIF 111 
COIF 111 
COIF 113 

coif 116 

COIF 115 
COIF 116 
COIF 117 
COIF 116 
COIF 119 
COIF 120 
COIF 121 
COE F 122 
COIF 123 
COE F 126 
C Of F 125 
COIF 126 
COIF 127 
COIF 126 
COIF 129 
COIF 130 
COIF 131 
COIF 132 
COFF 133 
COIF 136 
COIF 135 
COIF 136 
COIF 137 
CQEF 136 
COIF 139 
COIF 160 
COIF 161 
COE F 162 


SUBROUTINE INPTU Li ST) 

C This SUBROUTINE READS IN*UT oata 


INPT 1 
I NPT 2 
INPT 3 
INPT 6 


COUPON /CQNAF JhOB#N1»LOT#POUT#ITLEUO)#NAP|#DAY#NPG INPT 

1# JMUai*RRPN( lOltRCOUmm 111* INPUN#*OOPLOT#N#TRSN# BLADES INPT 

CJHNnNFCnN8/lRC0l#XRC0U3lf I RPP# C RP* ( 10 I # ZB AR ( 30 I # I YE 6 < 30 1 # SB t 30 I # INPT 

i ivEcnon$cnoitSNnmR6(3i)»RC(3iifTMnn f uT(30) inpt 

COPPPN/CQPC/VTPU 31>»EI6t30!»€lC!30>»G'H30)»THD<31> INPT 

COPPON /HINGIS/ LCH#LCHP1#LFH#LFMP1 # C HOF F # FHPFF # F CH# FFH IN*T 

1 #RPPPUN#COL*UN#LPM#i.PMPl#PHOFF#FPHf LOTS INP) 

COUPON/ COPJ/ XlP»PlP#ROPf POP# RBITA,XPSI»XT#XC INPT 

RE AL KIP»P!P#KQP#PQP#XB6TA»PP$I#PT#KC INPT 

LOGICAL l TVS# LOT# F OUT # TR$N# LOTS ► POOP LOT # IN PUN IN*T 

DIMENSION P ( 6 I IHP1 

NANE1-ST /INPUT/ NAPE#XT#POP#PtP#KlP#KOR#KBETA#KPSI# INPT 

1> C#N# JMUB#AZ8AR#TWIST#8LADIS#CHORO.CHOFF#FMOFF#PHOFF#PIAST# IRPP# INPT 
2RRPP# IRCOL#«COLL#Z#EI8.EICfWTPL#THO» EYE B# E Yl C • G J # SB# SC # RB # RC INPT 

DATA CVR/O. 0176533 / INPT 

C INPT 

C READ PROGRAP OPTIONS CARD INPT 

C INPT 

10 READ <5#90l> P INPT 

901 F0«PATU6#6X# A6*6X* A6»6X#A6»6X#A6»6X#A6) INPT 

IF (EOF ( 5 ) ) 320# 30 INPT 

30 CONTINUE INPT 

* • 0 INPT 


5 

6 

7 

8 
9 

10 

11 

12 

13 

16 

15 

16 
17 
16 

19 

20 
21 
22 
23 
26 

25 

26 
27 


38 
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965 


•Donor*. false. 

LOTS-. FALSE. 

LOT • .FALSE. 

•SOT • .FALSE. 

JNPUNa, FALSE. 

TRSN a .FALSE. 

LTWS a .FALSE. 

00 965 1*1.6 
IF(MII).EO.AHDECK) Ka 1 
IF(M(I).EO.AHNAAE) K*-l 
IFlM(I).Ea.AHMOOE) POUT*. TRUE. 

IF (H<I) .EO.AhALL”) LOTS*. TRUE. 
JFIM(l).EO.AHALLM) FOOT*. TRUE. 

IF (Fill .EQ. AHPLQT) LOT* .TRUE . 
IFIH(I).Ea.AHPUNC) INPUHa .TRUE . 
IF(HII).EO.AHTORS) TRSNa.TRUE. 
IF(MII).EO.AHNLTW) LTWS*. TRUE • 

IF (MU) .EO. AHSHAP) NOOPLOTa.TRUE 
CONTINUE 

IFIK.EO.O) 60 TO 319 


ORIGINAL PAGE IS 
OF POOR QUALITY 




IF K a o NO INPUT SELECTED 
K a 1 READ IN DECK 
K - -1 READ NAMELIST 

NPGaO 

IF (K.EO.I I 60 TO 20 
C READ CHANGES TO PREVIOUS CASE BT NAMELIST INPUT 
RE ADIS. INPUT) 

Nlatlal 
60 TO 55 

C READ INPUT DATA DECK 

20 CONTINUE 

READ! 5.902) N.'NE, ( ITLE II) » I-l. 10) 

902 F0RMATII6.AX.10if-) 

RE ADI5.90A) KT.NOP.MIP.KOP*KIP»KBETA»K»SI.KC 
90A FORMAT (8F 10 .0) 

RE AD 1 5. 905) N, JHUB. AZBAR. TWIST. BLADES. CHtlFF, F HOFF. PHOFF. PL AS T 
905 F0RHATI2I5.6F5.0. FlC.uJ 

RE AD 1 5. 2) IRPN. (RRPHI D.I-l.IRPH) 

RE ADI 5. 2) IRCOL.IRCOLLID.Ial. IRCOL) 

2 FORMAT I i5. 10F5.0) 

NJaH.l 

IFIAZBAR.E9.0.0) READI5.90A) 1 1 1 1 ). I«».N1 ) 


INPT 

26 

INPT 

29 

INPT 

30 

INPT 

31 

INPT 

32 

IMP i 

J3 

INPT 

34 

INPT 

35 

INPT 

36 

INPT 

37 

INPT 

38 

INPT 

39 

INPT 

40 

INPT 

41 

INPT 

42 

INPT 

43 

INPT 

44 

INPT 

45 

INPT 

46 

INPT 

47 

INPT 

48 

INPT 

49 

INPT 

50 

INPT 

51 

INPT 

52 

in* r 

53 

In;»t 

54 

INPT 

55 

INPT 

56 

INPT 

57 

INPT 

58 

INPT 

59 

INPT 

60 

iNP- 

61 

INPT 

62 

INPT 

63 

INPT 

64 

INPY 

65 

INPT 

66 

TNPT 

67 

INPT 

48 

INPT 

69 

INPT 

70 

. NPT 

71 


RE ADI5.90A) I E I d I D.I-l.N) 

RE AQI5.90A) (El It). I-l. N) 
READI5.90A) (WTPLI1). I-1.N1) 

IFILTWS) RE ADI 5.90A I (THDI I ) . I ■!» N1 ) 
READI5.90A) (ETEBI 1 1. I«1»N) 
READI5.90A) IETEC 1 1). 1*1. N) 

RE ADI5.90A) I 64(1). 1*1. N) 
READI5.90A) (SBII). t-l.N) 
REA0I5.90A) ISC ( I > a 1*1. N) 

RE ADI 5.90 A) (RBI I) .1‘l.Nl) 

READI5.90A) IkC I I ) > I • 1 » N 1 ) 

55 4HUB1»4HUB-U 
ZIDaO.O 

IFIAZBAR.EO.O.O) 60 TO 1A0 
DO 1AA Dl.N 
Z( 1P1)-Z(I)PAZBAR 
ZBARI I )«AZBAR 
1AA CONTINUE 
60 TO 1A5 
1A0 CONTINUE 

DO 1A7 1*1. N 
Z6AR(I>alim>-Z(I) 


IN i»T 72 
INPT 73 
INPT 7A 
INPT 75 
iNPT 76 
INPT 77 
INPT 78 
INPT 79 
INPT 80 
INPT 81 
INPT e2 
INPT 83 
INPT 8A 
INPT 85 
INPT 86 
INPT P? 
INPT 88 
INPT 89 
INPT 90 
INP1 91 
INPT 92 
INPT 93 
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147 

CONTINUE 

ton 

96 

145 

CONTINUE 

mot 

95 


IFUTMTI CO TO 160 

ton 

96 


00 190 I*1*JHUB1 

102 1 

97 

150 

THOI 11*0.0 

102T 

96 


THBATE«TMIST/2(N1 1 

loot 

90 


00 191 I*4MUB1,N 

ion 

too 

151 

TH0d*ll*TH0III*2BARIII*TNBATE 

102T 

101 

160 

CONTINUE 

ton 

102 


0*3 169 l* l*Nl 

I02T 

103 

165 

THIII*TMOdl*CVB 

tori 

104 


ifiinfuni cau cabos»","OOEno,nfs,b.oi 

ion 

103 


CALL STABT 

ton 

106 


BETUBN 

I02T 

107 

319 

CONTINUE 

I02T 

106 


HBITE(6,321I 

loot 

109 

321 

F0BNAT(6X,*7HNEITHEB DEC* OB NANELlST INFUT 0*T»ONS SELECTED 1 

ton 

110 

320 

CONTINUE 

IOOT 

111 


CALL CALFll 10.0,0.0,9991 

ton 

112 


STOF 

I02T 

113 


ENU 

ton 

114 


su about ine floutihasti flou 1 

HSU 2 

C THIS SUBROUTINE FBODUCES FAN HOTS • FlOU 3 

(••••••••••••••••••••••••••••••••••••• FlOU * 

COWON /COM/ J>«UB,NI,lOT,FOUT,ITLEdOt.NANC,OAV,NFC FlOU •> 

1, JMC c l.RRFNt 101 ,BCOLL I 3I»Z( 311 » l NFUN»*OOFL'IT»N»TBSN» BLADES FLOU 6 

CONH'IN/COHB/IBCOL ,XRCOL 131, I#FN»CBFNI l9t,IBABt)OI,EVE6(30),S9(10l,FlOU 7 

l ErEC(30l»SCI30lt$M31l»K6(31l,BCt3lt*TH<3lt,VT<3QI FLOU 6 

CONNON/CONO/CNATI9,9l»SONNATdOO, 31, IFLNI 100, 31 , INODE 1100,31* HOU 9 

1 NN3»N».9,CT( 31 l»ST (311, IB, 1ST, IBS 1 3, 10* 3l,tBFl3*10,3l HOU 10 

iOCICAL TRSN HOU 11 

OIHEMSION «»tO, VMM HOU 12 

ISIZE*! HOU lj 

CBU CALFinO.9,0.,-31 HOU 16 

XNI11*0.0 HOU IS 

XNI2I*B*FNI IBB" I *100, HOU 16 

XNAX ■ 6. HOU 17 

XNlN*.9*tFlXIIS.l-XRAXI*.9l HOU IB 


1X131*0,0 HOU 

IRFnC*RBFNd*FNI HOU 

DO 20 IT*1» 10 HOU 

I F H I V*200-IRFNC I .CT.OI CO TO 30 HOU 

20 CONTINUE HOU 

30 CONTINUE HOU 

IFtlV.EQ.3 .OB. 1T.E0.6 .OB. IT.EQ.9I tY*IT*l HOU 

tFiir.Eo.7i it*it*i hou 

XNI6I*IY*90.0 HOU 

YNUI*1.Q HOU 

YNI2I*FLAST *LOU 

NAXV*NIN0(e,|FIXIVn<2l*.01»0.9l-IFlItYNUI*. >1*0.913 HOU 

YNAX*HAXY HOU 

YM3I-0.0 HOU 

Yft<6l*6.*XNt6l HOU 

3 00 63$ 1*1.3 HOU 

IFIi.NE.3. ANO.BLAOES.NE .2.01 CO TO 639 HOU 

NBLi»B.60ES HOU 

call axesixhn, i.o,9o.o,yhai,yn(3i,yn<'.i,i.o,o.o, hou 

12HNATURAL FBEOUENCT, CBN, 0.19, 213 FLOU 

CALL CAlFlTlO. 0,1.29,-31 * FLOU 

IFII.EO.ll CALL N0TATE(3.7,0.10*TNAX, . 1?9,19HC0LLECTIVE NOOE.O. , 19FL0U 
II FLOU 


19 

20 
21 
22 
23 
26 
29 
26 
27 
2 » 

29 

30 

31 

32 

33 
36 
33 

36 

37 

38 

39 
60 
61 


40 



APPENDIX C 


5 


IMl.t'3.?) CALL F'0TATE(3.9.0.10«YMAV*.12S»UHCVCIIC NO0E*0.*UI 

IFII. E0.3t C All NCTATE(3«8*3.10 +Vnay* • 125* 13H5C ISSORS NO0E*O.*13l 
CET«3.A6-.258 RFlOAT ( IRCO l I 

CAU NUTATE (CET*R .5* .1* 17HR 30T COLLECTIVE • 17) 

CET»CET *1.657 
00 9 J-1.MC0L 

CAU nu**er(cet, 0.9. a, Rcoum* o.» it 

C£T«CET*.A29 

IFIJ. CE.lACeil 60 TO 5 

ORKHNAE PAGE B 
c.mrtNut OF POOR QUAMTV 

CALL NOTATE (CcT»6«5» • 1* 6HD£G«* 0*» 5 1 

IF(.NOT.TFSN) on TG 6 

CALL r«mm.36»8.65#3*I$m> 

CALL ROTATE ( 2.93* 0. 80*0.1* 7HT0RS I ON# 0.0* 71 
6 CONTINUE 

CALL FMmT<2.36*9.05*2*I$IZE) 

CAU NOTATE (2.69*9.0* .1* 1 IMHO* I Z *LANE»0.*11> 

CALL »itmm.36#9.2*»i»ism> 

CALL NOTATE (2. 76* 9*2# »1*10HV£*T *LANE#0.*10) 

CALL NOTATE (2.25*9.6#. 1*17hSTN NAK AN*l ITU0E#0. » 18 1 
CALL CAL*ir(2.25#J.375»3l 
CALL CAmrt2.*9,9.375»2» 

CALL CAL*LT(2. 69*9.375*31 
CALL CALUTO. 77*9.375*2) 

DO 0 

U«2.0MM*n*0.$8 

CALL NOrATE<XX#9.0#O.i*tTLE'(MI*O.O#6) 

6 CONTINUE 

CALL NOTATE (2. 39* 10*00* • 15# 33HC0UFLED BLADE NATURAL MEOUENC I E$r 0 . 
1*33) 

VNANE*NANE 

CALL NOTATE (0.5*10.0»Q. 1# 6HC AS E* 0.0*6) 

CALL NUNBE9I 1.0* 1C. 0*0. l#XNA«E *0.0*0) 

CALL NOTATE « 7*63* 10 .00* • i* 10H( )*0.*10) 

CALL NOTATE! 7*656# 10*00# . 1#DAY *0.*9| 

CALL CALNLTCXNIN*-0.29»*3) 

CALL AXE$(.0#*0#.0#XNAX»XNt3l#XN(6)#l.*.0*i5H*OTO* SNEED#** N* . 35* 
1*151 

00 200 IFF*1#8 
VSAQT-IFF 

IECI.NE.il GO TO 210 
THIS BATH FOR COLLECTIVE NODES 


62 

63 

66 

65 

66 

67 

68 
69 

50 

51 

52 

53 
56 

55 

56 

57 

58 

59 

60 
61 
62 
63 
66 

65 

66 

67 

68 

69 

70 

71 

72 

73 
76 

75 

76 

77 

78 

79 

80 
81 
82 
83 
86 
85 




I F t*ODt IFF#NBLD)»NE«0) GO TO 210 

RLOU 

86 



GO TO 236 

RLOU 

87 


210 

IF ( I.NE.2) GO TO 220 

RlOU 

88 

C 

TMI« RATH for cyclic .does 

RLOU 

89 



• F (NODI IFF.H81DI.E0. 01 CO TO 220 

RlOU 

90 



CO TO 236 

RLOU 

91 


220 

Ifd.NE.3l CO TO 200 

RLOU 

92 

c 

THIS RATH FOR SCISSOR NODES 

RLOU 

93 


2 36 

CONTINUE 

RLOU 

96 

c 



RLOU 

95 

c 


RIOT 1* "2* 3 RER REV LINES ON RLOT 

RLOU 

96 

c 



RLOU 

97 



XRO$*XNAX 

RLOU 

98 



YROS-YSROT 

RlOU 

99 



If (YSROT.lt. YHAXt CO TO 233 

RlOU 

100 



VROS*YNAX 

RLOU 

101 



XROS>YRO$/(0.23*YSROn 

RlOU 

102 


2 35 

CONTINUE 

RLOU 

103 



CALL DA$HLN( 0.0* 0.0* XROS* YROS *0.21 

RLOU 

106 



CAU NUNBER( XROS* VROS-O. 05*0. 1* YSR0T,0.0»-1 I 

RLOU 

105 



NOTATE (XROS* 0. 08 37* YROS-O. 05*0.1* AH/REV* 0.0* A) 

RlOU 

106 



CALL C ALRLT(0* 0* 3 I 

RLOU 

107 


41 
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200 CONTINUE 

C ALlT ROTO* NATURAL FREQUENCY 

00 410 l9*l»tRRN 

oo aid *iST»i*hccl 
J1*!4$M'T» IS»U 
J2* Iif < 1ST* 19* 1 1 
IMJ2.tT.Jl> GJ TO 410 
00 AGO JOl.J? 

«»i*lno»ii 

*•*«»«< 

y*$onnatu» nsru 4»*y*m 
call RNmru»Y»**isuc> 

AGO CONTINUE 
A 10 CONTINUE 

CALL C4L*LTI9.0~MNIN*~1.0*~31 
435 CONTINUE 

CALL NFRAHg 
*6 TU«N 
END 


non ioa 

A LOU 100 

nou no 
nou ill 
nou ik 
nou m 

A LOU 114 

nou as 
nou no 
nou nr 
ncu no 
nou 119 

8L0U 170 

nou Hi 

nOU 172 
nOU 123 

nou 174 

nOU 123 
nOU 126 


c 

c 

c 

c 


SUBROUTINE ST At T STAt 

• SfA« 

This 5U9AOUT1NI MINTS our IN»Uf OATA A 90 RECASTS data TO STAR 

LU»nO MSS AT STATION 80UNOARY RE*RE$f NT AT IGN STAR 

•♦♦♦••♦••••••♦♦•♦♦•••♦♦♦•♦♦♦♦••♦♦♦♦♦♦♦♦♦•♦♦♦♦♦♦♦♦•♦♦♦•••♦A STAR 

CGN*TN /COM/ JHI 9»Nl»LOT»ROUT»lTL£UOUNANr,0AT»NFG STAR 

l»JHUSl»«tRNI10ltRCOU(3l»2( )ll>lN7L«N v NO0RLOTtN.TtSN»iLA0€$ STAR 

C^NN3n/cnM/lRClL#ARCQLl3>»tRRN#CRpN(10l*7RAft(30>*iTCSf 30 > * SB < 30 > *$T A* 

i e tec 1 30 1 rc<so»»SNmi»Rsnn*itcm>»TMni»»uT(30i star 

CG*«ON/ CbNC / HTRl(3ll»Elft(30>»£ICf90l»£J<30 I *THQt 31 1 STAR 

C3NN0N /HINGES/ ICh,lCHR 1#LFM»IMR1 *CMnfF* FM0FF *FCM*F FH STAR 

1 . «MPU*,C Utl'N.l fM.LtNtl , RnOf F,T/m. LO*S STAR 

CONNON# ONJ/ MnNtt,R0nN3».*S£TA*Rm#RT*Rt star 

C >n*ON/CONT/ ET9t31)»EYCf31l*EYft(3l)*f av RM31 I# EMC (31 T#CNRRI 311* STAR 
1 €Ntt{31l»€NRT( 31 )»6NBSWl31l»ENfRtfl3l 9»THHOni)»ENftRWl31) STAR 

real * l A# N l R# k0R# HOP* ABET A# a R$ l **T**C STAR 

LOGICAL TR$M STAR 

DINENStON C F I 301 STAR 

DATA CVR*C¥«R$/ 0.0174533*0.1047199/ STAR 

R*Z<NII STAR 


1 

2 

3 

4 
3 
6 
7 

a 

9 

10 

11 

12 

13 

14 

15 

16 
17 
ia 

19 

20 



ranaun-RRRni IRA*! 


STAR 

21 


IM IRRN.CT.7I RtRRUfc*RRRftt IRAN-1 1 


STAR 

22 


CJl RUN* PC PL l ( I PC 01 1 


STAR 

23 


Iff IRCGL.GT.21 C0lRUN««C0l L C IRCOl-1 1 


STAR 

24 

c 

CONFUTE 9LA0E ANGLE AT ENO Of EACH SFGNfNT 

KCOU IN DE6KECS 

STAR 

25 

c 

XCOl IN RADIANS 


star 

26 


03 7 0 1*1* I RCGL 


STAR 

27 

710 

vRcoim«RCom n *ct« 


star 

28 


00 750 1*1*1 RAN 


STAR 

29 

750 

CPRNCI 1-PRAFM IMCVPAS 


STAR 

30 

C 

CONFUTE BLADE STATION f OR r L A At JNG MING* * 

CMOKONISE HIN6E. AND 

STAR 

31 

C 

A ITCH HORN Cf f SET 


STAR 

32 


*UN7*0.0 


STAR 

33 


DO 974 I ■ 1 * N 


STAR 

34 


K*l 


STAR 

33 


$uhz***jhz»:bar! n 


STAR 

36 


IFfSUNZ.GT.RHGFF 1 GO TO 973 


STAR 

37 

9 74 

CONTINUE 


STAR 

39 

875 

CONTINUE 


STAR 

39 


LPHA1-K 


STAR 

40 


L RM«K-I 


STAR 

41 


F»H* ( PHOf f*SUNZ«7BAR( L AHA 1 )) /ZB AR ILAHAI) 


STAR 

42 


5UNZ*0,0 


STAR 

43 


42 




00 876 t»UM 


star 

99 


SUN{*SUN2»|0M<tl ® 


STa* 

STAR 

95 

96 


ifisuni.gt.fmoff) go to «77 ‘ 'P rUOR QUALffr 


STAR 

97 

ih 

continue 


STAR 

98 

877 

CONTINUE 


STAR 

99 


IFMP1N 


STAR 

50 


IFNN*1 


STAR 

51 


FFINIFNOFF'SINMIAIIIFWUI/IIMUFIPII 


STAR 

52 


SUN2*0.0 


STAR 

53 


00 OTO 1*1, N 


STAR 

59 


■ •I 


STAR 

55 


SUN2*SUNZ»ZBMIl> 


star 

56 


IFftUNZ.GT.CHOFFI CO TO 079 
CONTINUE 


STAR 

57 

in 


STAR 

58 

879 

CONTINUE 


star 

59 


ICNP1>* 


star 

60 

jf 

ICH*R-1 


STAR 

61 

m 

FCH*(CNOFF>SUNZ»ZBAR(UHPin/ZBARILCMPll 


STAR 

62 

* c 

CONPUTE CENTRIFUGAL FORCE, BLADE NA*S.ANO BLADE FLAPPING 

INERTIA 

STAR 

63 


FLPINT-HTPL!N1>*(«-FhOFFI**? 


STAR 

69 


TOTNAS*HTPL INI | 


STAR 

65 


00 09 1*1, N 


star 

66 


CFt 11*0.0 


STAR 

67 


00 07 IJ-I,N 


STAR 

68 


CF 1 II*CF I I I »HTPl 1 IJ)RZ4AR(tJl*IZ(IJt*Z(IJ*I>l*.9 


STAR 

69 

•7 

CONTINUE 


STAR 

70 


CFni*ICFIII»MTPl(Nll*«l«CVRPS*C0RPS/306.A 


STAR 

71 


IFII.lE.LFHI GO TO 09 


STAR 

72 


TOTNA$-TOTNAS*MTPl(t!«ZBAR(tl 


STAR 

73 


FLPINT-FLPINT**.*! PllI»*ZBAtIII*IZ«I»*ZCI,l»-2.*FMOFFl**2»/C.O 

STAR 

79 

85 

CONTINUE 


STAR 

75 


FI PINT *FL PINT/I 32 .2 9166. 1 


STAR 

76 

C 7RIHT SKIT IMUT • 


STAR 

77 


NPG*NPG*1 


STAR 

76 


HR ITEIb»909) NPG, DAT, NAME, (ITLE(A),J*1,10> 


STAR 

79 

905 

F0BN*T<1H1,18X,6HPAGE,I3,2SX,36HNV«LF'T60 ANALYSIS ASA 

T* 78670 

STAR 

80 


1 ,29X, A10,//10X,SHCA$E , Ib,23X, 10A6/1 


STAR 

81 


HR1TE(«,906| 


STAR 

82 


DO 906 1*1, N 


STAR 

83 


HR ITE 16,9101 I,ZBAR| I l,EIB( 1 1, EIC (I ),HT*L ( I l,TH0 1 1 1, CFlII 

» E YE 8 ( 1 1 » STAR 

89 


1ETECIII,RB(II,RC(II 


STAR 

85 

910 

FORNAT (I6,FB.2»2E11.6,F8.6,F8.2,EI1.6,tX, 2E13.6, 29X,2F8.3) 

STAR 

86 


IF ( .NOT .TRSNI GO TO 906 


STAR 

67 


HR ITE 16,908 1 Gjm,SBII),SC(II 


STAR 

86 

908 

FORNAT ! 1H»*87X» E 13.6, 2F8. 3) 


STAR 

89 

909 

CONTINUE 


STAR 

90 


HR ITE (6, 9091 HTPL(N1),ThO(N1I,RB(N1I,RC(N1) 


STAR 

91 

909 

F0RNaTI3X,3hTIP,27X,F9.6,FB.2,67X,2F8.3/> 


STAR 

92 

906 

FORNAT I 6X, 6HSEGNNT, 7X, 8H El , 9X. 9HHT/ IN, 6X» 5HTHIST, 5 X 

#9HC«F •» 

STAR 

93 


1 10X, 3H I BB, 10X, 3HICC,8X,BH GJ , 3T, 12H6MF AR OFFSET, 61, 

11MC.0. 

ostar 

99 


2FFSET/7X, 6HLNTH,6X, 6HBE AN,9X,9HCH0R0, I3X, 6H( 1NB0 l,6X,6Hl IN60),66X 

star 

95 


3, 6HBE AH,6X, 9HCH0R0, 3X,6H8E AN, 3X, 9HCH0R0/7X,6H IN,6X,16H* 

- i6-kn**star 

96 


92 —f5X#5HtB/IN#5I#3M0ECO8f 9HLB/8786#2 > 9»ol5HIN-lB-SEC9*2/IKf 8X 

• star 

97 


9BHLB-IN**2,3X,2HIN,6X,2HIN,6X,2HIN,6X,2HIN) 


STAR 

96 


HR I TE( 6, 882 » BLADES,R,TOTNAS,JHUB,FLPINT 


STAR 

99 

862 

FORKATt 10X»F6.0»7H BL AOE S# 181* 6HR AO IUSo FB . 2» 3« 25X# 12HBI AOE WE 1ST AR 

100 


1GHT,F6.2,3H LB/10X, 18,13H HUB SEGHENTS,96X,18HBL AOE FLAP 

INERTIA 

#$TAR 

101 


2 F 8.2, 1 1H $LUG-FT**2/> 


STAR 

102 


HR ITE (6, 883 ) FH0FF,K6ET*»H0»»*QP,Ch0FF,KPSI»HIP»IUP»PH0FF 

,RC,KT 

STAR 

103 

863 

FORNAT! IX, 13HFLAP HNG OFST,F8.2,3H IN.9X,11HFLAP SPRING, E9. 3, 10H 

ISTAR 

10 A 


1N-L8/RA0,5X, 16HHUB O.P. INERTI A, F 6. 3, 1 3H LB-SEC**2/IN, SX, 


STAR 

105 


2 16HHUB O.P. STIFF, E9.3.6H LB/ IN/IX,! 3HLLG HNG OFST ,F8,2,3H IN, STAR 106 
3$X, 11HL AG SPRING »E9.3»10H IN-LB/IU0»5X,16HHUB I.P. INERT I A, F6 . 3, 1ST AR 107 
63H LB-SEC**2/IN»!>X» 16HHUB I.P. STIFF, E9.3,6H LB/IN/1X, 15MPI TCH HORSTAR 108 
9N OFST , F6.2, 3H IN.5X. 11HCNTRL STIFF, E9. 3, 10H IN-LB/RA0/26X, STAR 109 



r> o n 


WRMHX C 


ft ITMNAST TMSlTt JN *E9«3*10H IN-LB/RAO/I 
C CALCULATE COEFFICIENTS INDEPENDENT OF COLLECTIVE ANGLE AND ROTOR I 
0) SOS I«2*N 

SEGIN-O.SFUTPvI t-l»«2MR( I-l l/3*6.A 
SCGQUT«O.S*UTPl(ll»2ftAR(l)/)8ft.A 
SN(ll»SE6lN*$EG0UT 
(Mil ll«SEGlN*R8tt-ll*$S60UTFRBl II 
ENRC(tl*SE6lN»RC(t-ll»SE60UT«RC(I» 
Enpr(1I»SEG1n*RB(I>U««CII-1IaSE60UT*RS(IIFRC(1I 
Eriui*t.)«i(i(«i!*u>tMMMi»mMn»mMm«i((iMiiiMii 
l»F2/l2.0*Re(I-ll*»2lASE60UT*12BAR(U**2/l?.0*R8(II*F2l 
CfCtII*0.3*(EVECII-ll*2BAR(t-mEV£C(t>FZftAR<III«SE6INF(ZftAR(l 

iMi/u.o*icii*UM;i»s((aurMtu»in*M/iMi»(cui«*ii 

SOS EYRlI»»EYi(tl*EYClII-ISEGINFTftAR(l-l»**2*SE60UTFt9ARIII**21/6, 
S<* (1 1*0. 3»UTPLlll*2BA.t( 11/386. A 
SN(N1I*Q.S*HTFI(NI*26AR(NI/3B6.A 

EMIllliSRIllMIIU 

ENRB(Nl l*SN(Nt I 9RB(NI *HTPL (N1I*RB(N1I/306»A 

£nrc(ii*shiii*rciii 

ENRC(N1I*SN(N1»««C(NI*NTPI(N1I*RC(N1)/3B6.A 

ENRRtll«SN(l)*ACtll*RB(ll 

ENRR(N11*SN(NII«RC(NI*R0INI«UTPUN1I*RC<N1I»*B(N1I/386.A 
EVB(1I«0.3»(FyEB( 1I*2SAR(1|| ♦SN(llF(rNAMll**2/12.0*R8lll**2l 
£rG(ftll*0.5F(EYE8(ft|A2BAA(Nll*SH(NHMRB(NlFF2*2ftARIN|FF2/12.C 
l ♦HTPl(Nll*RB(NH**2/18ft.A 

EYC(lt*6.3*(EYEC(l)F26AR(ll l«SN(ll*(ZBAR( U**2/12.0»«CU>»*2> 

E YC(N1I*0.SMEYEC (N 1*18 AR(NII*SN(N1>* IRC (N|*F7*2BAR(N»F*2/12.« 

inikitMcmiiMmit.i 

EYR(ll*EY8(H*EVC(ll-SNIl|F2RARUI**2/6.0 

FVR(Nll-FVB(Nll*EVCtNll-SN(Nll'2BAR(N|FF?/6.0 

SN(Nll*SH(Nll«HTPl(«ll/3R6.A 

RETURN 

ENO 


SUBROUTINE SUNNY SUM 1 

SUNN z 

THIS SUBROUTINE FRINTS OUT A SUNNARY OF NATihaL FREQUENCIES • SUNN S 

ft************************************************************ SUNN A 

CONNON /CONA/ |HUB*N1 *i.0T»P0UT* I TIE 1 10I*NANE»0AY»NPG SUNN S 

1» JNUBl»RRPNI10l»RCOlLCII*21311.lNPUN*NnO*LOT»N,TRSN*BLADES SUNN 6 


C0NN0N/C0NB/ IRCOltXRCOl (Sir IRAN, CRFNt 101* 2SAR< 30I*E YEB( 39I»SB( SOI* SUNN 


I t TEC (30t*SC(30l»SH(31l»RBI31l»RC(31l*TH(31l*HT(3CI SUNN 

CONNLN/t OHO /CHAT! 3*5 l*SONNAT( 100* 31* 1 Ft. N( 100*31* INODE ( *00* 3 ** SUNN 

1 HN3*HN3»CT( 311 *$T( 31 1 * IB* 1ST* IBS ( 3* 10* 31 * IBEI 2* 10* 31 SUNN 

CONNON /HINGES/ LCH»LCNP1»IFH,LFHPI *CHOFF*FHOFF*FCH*FFH SUNN 

1 »APHPUN,COLPUN»LPH,LPNPl»PHOFF,FrH» LOTS SUNN 

0 I NENSION FINE ( A* l SUNN 

DATA *INE /6HVERT >6HH0RII *6HT0RSN *6H / SUNN 

NP6*NP6*1 SUNN 

UR ITE (6*001 1 NP6.0AY*NANE»(ITLEUI*J-1*10I SUNN 

901 F0RNATI1M1»18X»AHPAGE»I3.23X.3AHHYKIESTA0 ANALYSIS NASA TN ?**Tn SUNN 

1»23X,A10/1SX.SHCASE » Ift* 1BX* 10A6/ I SUNN 

IFIBIADES.6T.2.0I 60 TO 10 SUNN 

UR ITE 16*9091 SUNN 

90A FORNaT ( 1 3X» 19HC OllECTIVE* ?TV» 13HS C I S $ 0 R S * SAX* SUNN 

1 11HC Y C L I C .// SUNN 

1 6X« 31 AOHNATURAL ROOT ROTOR NAX NUNBER * 3X1 » /TX, 3 (SUNN 

133HFREQ COLL SFEEO INERTIA OF * «X I » /7X» 31 37H/REV 0E6 SUNN 

2 RFH PLANE NODES *6X1 I SUNN 

GO TO 13 SUNN 

10 CONTINUE SUNN 

HR ITE (6 *903 I SUNN 

903 FORNAT I 59X* 13HS CISSORS*/ / SUNN 


7 

B 

9 

10 

11 

12 

13 

1A 

15 

16 
17 
IB 

19 

20 
21 
22 
23 
2A 
23 
26 

27 

28 
29 


44 



APPENDIX C 



X *■»*. tOMNtTURtl toot ROTOR Nt« NUNRER t 

SOT* SUNN 

30 


UDNEREO COll SREEO INERT It Of /*0», STH/RE V 

DE6SUNN 

n 


» Rtw RltNE NODES /> 

SUNN 

32 

l* 

CONTINUE 

SUNN 

3) 


LINES-0 

SUNN 

3* 


lNt-1 

SUNN 

3* 


1N2*1 

SUNN 

>0 


|NJ»l 

SUNN 

3? 


DO 60 IST«X» IRCOl 

SUNN 

30 


lt*E$«t INES*l 

SUNN 

10 


NRITf <6.4021 

SUNN 

*0 

902 

FQRNtT (1H 1 

SUNN 



SO so 18« l » IRtN 

SUNN 

M2 


NO-INEtlET. 1 <*.*>- IBS! 1ST. 18. J» *X 

SUNN 

Ml 


IFmtDES.ST.2.D> GO TO TO ^HIGINAL PAGE 

SUNN 

SUNN 

MM 

MM 


«8.i#£UST.is.n-i»s«iST.!6,n.i OF POOR OITAirrv 

NC«IBEIIST»18»?l-leSII'T.IB,2)»X MUAAJll 

SUNN 

SUNN 

66 

Ml 


Nt«NtK0IH8.NC.NDI 

SUNN 

68 

20 

CONTINUE 

SUNN 

69 


IMNt.EO.OI GC TO SO 

SUNN 

50 


tMl!NES*Nt.iT.*OI GO TO 21 

SUNN 

51 


NtG»Nt6*l 

SUNN 

M2 


W* ITT 16.601 1 NtG.DtV.NtNE.ITlE 

SUNN 

53 


IF<8ltOES.GT.2.0> GO TO TX 

SUNN 

55 


WRITE (6*906 1 

SUNN 

55 


GO TO 2? 

SUNN 

56 

21 

CONTINUE 

SUNN 

5? 


WR ITE 16. 60S I 

SUNN 

55 

22 

CONTINUE 

SUNN 

50 


IINES-O 

SUNN 

60 

23 

CONTINUE 

SUNN 

61 


RRITEI6.60TI 

SUNN 

62 


IINFS«IINES*X 

SUNN 

63 


• WtXOINft.NC.NOI 

SUNN 

66 


l INES-l INES*"t 

SUNN 

65 


DO 60 I • I . Nt 

SUNN 

66 


WR ITE 16.6021 

SUNN 

6? 


IMI.GT.HDI GO TO 2S 

SUNN 

68 


ENtT«SONNtTl INI, 31 

SUNN 

60 


IF(RRtNUB).NE.O.O) FNtT»FNtT/RRNNCIN» 

SUNN 

70 


it«iHNiiN).ji 

SUNN 

n 


WRITE 16.6061 FNtT»RCOll(ISTI»RRtN|ItT.tlNECltl.I NOOE (IN}.)) 

SUNN 

72 

90e> 

EORNtTI IM..19X.F10. I.TFT.l. 11.66. Ill 

SUNN 

73 



I N 3* I N 3* 1 

SUNN 

7* 

25 

IMHA0E$.GT.2.0> CO TO 60 

sunn 

75 


If < l.GT.nil CO TO 10 

sunn 

76 


f NlUSOWNATlINUl ) 

Sunn 

77 


X H*«M< 16) .*€ *0.01 fNAT*SO**NATUNltl >/RR6nt tRI 

sunn 

76 


I’MPIMINWI) 

Sunn 

79 


WRITE (6»907| fNAT»RCOU (IST)f RRM< lBt*nNCt!6)»lN00!tlMl*n 

Sunn 

60 

007 

f ORBIT UH4i '*«F10.1»2f7,l»lY»A6»n> 

sunn 

61 


WOMN1M 

sunn 

62 

10 

tfU.GT.NO GO TO 60 

sunn 

63 


fnat-sonnatun:*: 1 

sunn 

66 


If <RRMU8> .HE .0.0) f NAT«SO*N*T tl N2 » 2 1 /RRM U 6 > 

Sunn 

65 


Uh2»2J 

Sunn 

66 


WRITE (6.906) fNAT.RCOUUST) .RRMUBUnnHlMtlNOOH !N2»2> 

Sunn 

67 

906 

f ORfUTl lM*.88**f 10. )#2f 7. W 11# 66 » lit 

Sunn 

66 



sunn 

69 

60 

CONTINUE 

sunn 

90 

50 

CONTtNUE 

sunn 

91 

60 

CONTINUE 

Sunn 

92 


RE TURN 

sunn 

93 


END 

sunn 

96 


45 



APPENDIX C 



SUBROUTINE NOOSHAP|9,COU-R»N,FREORR.N,«,«ENO> 

WOOS 

1 



WOOS 

2 


THIS SUBROUTINE PLOTS the node shape 

WOOS 

3 


CJNHJN /COHA/ JHtB-Hl-LCT-R JUT • I TIE ! 10 1 -NAHE- 0AT-NP6 

NODS 

6 


1- JNU81-RRMI 10)-RC0Lli3).H3ll» INPON-HOORlOT-N-TtSN- HADES 

WOOS 

5 


LOGICAL T»SN 

noos 

6 


3 1 HENS ION BI3-31I-RI T3> *6£ AH03 T*CMO»OT J3 >»TO*SI ONI 33 )» NODE IS) 

woos 

T 


DATA H'JOE/lOHCOLLECTIVE.iOMCTCUC . lOHSCISSOtS / 

WOOS 

0 


XJRIG-0.0 

WOOS 

9 


XSCALE*0.2 

WOOS 

10 


AIENTH-1.0/1SCALE 

WOOS 

11 


YORIG*-l.O 

WOOS 

12 


YSCALE*0.50 

WOOS 

13 


30 10 I *1 • HI 

WOOS 

19 


R(II*2<It/2<Nl) 

WOOS 

15 


8E AN| I )*8 ( l - 1 ) 

WOOS 

16 


CmORDII)*8I2.I) 

WOOS 

IT 


TORSION! I )*3I3, I) 

WOOS 

10 

10 

CONTINUE 

WOOS 

19 


R INI* l) *XQR 16 

WOOS 

20 


><N1«?I*XSCALE 

WOOS 

21 


8EAH!N1-1)*Y0RIG 

WOOS 

22 


9EAR|Nl*2)*VSCALE 

WOOS 

23 


CHORD!Ni*l)*TO»IG 

WOOS 

29 


CH3P3(N1«2I-VSCALE 

WOOS 

25 


TORS I0NIH1*! 1*TURIG 

WOOS 

26 


T0»SI0N|Nl*2)*rSCAlE 

WOOS 

2? 


t*N0DI*»2) 

woos 

26 


IFII.GT.O) CALL CALPlTI 1.0-6. 5.-3) 

WOOS 

29 


CALL IINPLT(R»9EAH-n)*1»1-1»1*0) 

WOOS 

30 


CALL IINPlTIR»CHORO-N1- 1-1- 2-1-0) 

WOOS 

31 


H-0.15 

WOOS 

32 


IFIRPN.GT.O.O) gc TO IS 

woos 

33 


CALL NUHBERf 0. *»*.>» M,FRt0PR»0. 0» 01 

WOOS 

39 


CALL NTTATEI1.7-A.9-H-AH CPN-O.O-A) 

woos 

35 


GO TO 16 

woos 

36 

15 

CONTINUE 

woos 

37 


CALL NUN8FR 10. 8- A. S, H, FREQRR- 0. 0- 3 1 

woos 

30 


CALL NOTATE! 1.6-a.6,h,6.< / REV-0.0-6) 

WOOS 

39 

16 

CONTINUE 

woos 

90 


CALL NOTATEI2.S-A.3-H-NOOEIN)-O.0-I0> 

woos 

91 


»scale*xscale*i-i.o: 

woos 

92 


CALL NUH8FRIO.9-A.3-H-RRH-0.0-0) 

woos 

93 


CALL NOTATE! 1.6-A.3-H-3HRRH-0.0-3) 

woos 

99 


CALL HUP8ERI0.8-A.l-H-C0lL-0.0-l) 

woos 

95 


CALL NOTATE ll.6-A.l-H.AHCOLL-0.0-A) 

WQDS 

96 


CALL AXES IUENTH, 0.0, 180.0. XLENTH. 1.0, XSC ALE. -1.0- 0.0- 6HR AO IUS-H-6HOOS 

97 


1 ) 

woos 

90 


CALL AXES 10. 0-0. 0-90.0- A. 0--1.00-YSC ALE--1. 0-0.0# 10H0EFLECTI0N-H- 

woos 

99 


110 ) 

woos 

50 


CALL PNTRITCA.2-A.6.1-1) 

woos 

51 


CALL NOTATE lA.A-A.S-H-AHVERT-0.0-A) 

woos 

52 


CALL PNTPLTIA.2-A.A-2-1 ) 

woos 

53 


CALL NOTaTEIA.A-A .3-h-SHHOR IZ-0.0,5) 

woos 

59 


IFI.NOT.TRSN) 60 TO 20 

woos 

55 


TSCALE*YSCAIE*«-1.0) 

woos 

56 


CALL L INPLMR, TORS ION, Nl, 1-1-3, 1,0) 

woos 

57 


CALL RNTPLTIA.2-A. 2-3-1) 

woos 

50 


CALL NOTATE! A. a-a.i-h-7hT0RSI0H-0.0-T) 

woos 

59 

20 

CONTINUE 

woos 

60 


IFII.GT.O) CALL CALPLTI0.0--S.S--3) 

woos 

61 


IFIX.EO.A END .OR. I.EO.O) CALI NFRAHE 

woos 

62 


RETURN 

woos 

63 


END 

woos 

69 
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APPENDIX D 


ORIGINAL PA^ 
OF POOR Qt 4 


LIBRARY SUBROUTINE MATINV LISTING 


The Htting for the library subroutine MATIKV follows! 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C 

c 

c 

c 

c 

c 


XUfckOUT | ML HA 1 t Nv tMAX ♦ N% A % M % B * 1 OP * 0 t t LNH ♦ IhCA^t « IP I VOt » IKM 


PU^POSt - HAT INV lNVC«Ts A HLAl SvXjAk't HATWIX A* 

IN At'OmON THt MOO TlNk. SOLVt S THt HAIP|X 
CQOAMON ft IS A MAT W I X OT CONSTANT 

vactowsb r*n pt is ai so an option to mavt twr 

or Tt'WH INANT tVAlUAUO* IT 1«t lNVr«sr IS NOT 
NtiOtO* UST OU. IH TO SOt Vt. A SVSTFM 0* S I Nt »t T ANt OPS 
feOUATlONS A NO iHfTTAC TO r, V*luAT» A UTTlRMlNANl 
► Os' SAvlN^ M M* A NO s*OWAx»t • 

- v Al t. HAT INV lMAX*N% A*M«tW ! OPtlH 1 1 PH * I SC At t * IP, VOT « l*K t 

HAX - THt MAXIMUM ONOtS Of A A \ j» T A T t D IN THt 

0 I Ht NS I ON STATtMtNT 0* tML LAcc I NC» »*NO v *«AM« 

N - T^t Oppt <4 Of A % I « i t • N • L t • M A X • 

A -A T*0-0|TH NolONAL AHWA V Of • Mt CCH f T K 1 1 NTS* 

ON Nt U*S\ 10 TM* CALLING PPOUWAM* A iNVtNSt 
IS O IN A* 

A MSvl \U OlHi.NsI ONt.0 IN THt cALclNv* ONOOHAM 
*Mo Mh'sT OlMtNMON H Ax A No SI i ONO P I Ml Ns 1 ON 
AT It AST N* 

M -THt NOHOtK O* xOtwMN VCCU)NS IN ti« 

**v> SlONAi THAT THt SoPSSXM INt IS 

w»M V sou I t > ► ON' INVi s's | PNtHOfct vt s« 

|N I ML v Alt St AT t M| M AN t NtKY v OK'K't - 

SPONOI Su Tv' o MSI Pf PNt St NT* 

ti - A TlfcO *P I M| N,- 4 I ONAu ASK A > O* THt vv'NsTANI 

yt v 10s' 0% ON N? IpK'N T O V At l. I NU PKOhNAM* 

X I*' Mos't l' IN O* P smiH»i O oAvt Ms ► 1 WS 1 

OMUNSlON MAX ANO I Is St v ON* AT i I AM M • 


lOP vOMPutl Ot UHHlNANl OPTION* 

100*0 vOMI'olts TNI HAimx iNVtWSI ANP 

Of Tt NH I NANT • 

lOO»t vOMOv»Hs f Ml MAI Six INVI PSt v'Nl X • 

01 U s'H ►ON I0O*v-IN w'NJt»Nx T I ON fcllH ISCAU 

NtONtstNIS THt vAltU Of THT MTtN'MlNANT 
C C'T A, i't ’|AlrV. S’Uv'O'. 

V, Ot U A > UH I * K'M | i I 0 • * l OO i l Sv Ai l n 

c t tit v O v OO 1 A I l ON tU I lA I smOKS-P Nv' I Ot 

C At Tt Mono IN T n\ osf k‘ »S\'i.KAW s | Nv t I* 

^ 1 M( 0*V0l K v'> A | '* i Ah’vA ►% ANi' v'K I »4l 

v MAi.MToOt v'T MS III MINI'. AM t AKv<| i SMAl l \% 

t 1 Mt lM M A i v Al x i v A T l ON M A > v Ai 1 0 i'\ » s'* i i' * 

v. * v <NOt t O* ' • 01 1 ■ ‘"'M '»{ 1 tv'' » t ^‘v x * OO 

v SlNvl.A.V VAh. |\ vONOMK'N* * 0*1 ' l T K' 

v |vO I , O'i v' • SMv'vi 0 ot v*uvV, o OX OSOvN'AMt S' 

ON M IvS'N U" MAIN PKvVvAM, 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

C 


c 

c 

c 


I scale - a scale facto* computed by the 

SUBROUTINE TO AVOID OVERFLOW OR 
UNDERFLOW in THE COMPUTATION OF 
THE QUANT I TV *DC TERM* 

IPIVOT - A ONE DIMENSIONAL INTEGER ARRAY 
USED BY THE SUBPROGRAM TO S'ORE 

pivotol information* it should et 

DIMENSIONED AT LEAST N. IN GENERAL 
THE USER DOES NOT NEED TO MAKE USE 
OF THIS ARRAY* 

IWK - A TWO-DIMENSIONAL INTEGER ARRAY OF 

TEMPORARY STORAGE USED BY THE ROUTINE* 

IWK SHOULD HAVE ITS FIRST DIMENSION 
MAX* ANO ITS SECONO 2* 

REQUIRED ROUTINES- 

REFERENCE -fox*l* an INTRQOuC T 1 ON TO NUMERICAL 

linear algebra 

STORAGE - 542 OCTAL LOCATIONS 

language -fortran 

fortran 

library FUNCTIONS -abs 

LATEST REVISION - JULY 1973 -CMP6 

MATRIX INVERSION WITH ACCOMPANYING SOLUTION OF LINEAR EQUATIONS 

DIMENSION IP I VOT <N> * A <MAx*N ) (MAX *N ) * I WK (MAX. 2 > 

EQUIVALENCE UROW'JROWI' < ICOLUM* JCOLUM) • < AMA X % T. SWAP ) 

INITIALISATION 

!SCALE*0 
R1 =10*0*«100 
R2 = 1 • O/W I 
DETERMa I *0 
DO 20 I * N 
20 IPIVOTI J)=0 
DO 550 I = I *N 

SEARCH FOR PIVOT ELEMENT 

AM AX * 0*0 

00 105 J*l «N 


IF 

< IPIVGTC J)~l 

) 60 • 

105« 

60 

DO 

100 K-l |N 




IF 

( IP I VO l (N >- l 

) BO* 

100' 

740 
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00 IF <ABS(AMAX >-AdM A <J*K ) ) »85*lOQ* 100 
85 IR0W*J 
JC0LUM*K 
AMAX*A< J,K ) 

100 CONTINUE 

105 CONTINUE 
IF (AMAX) 2 10*1^6010 

106 OETERM*0*O 
1 SCALE *0 
GO TO 740 

1 10 I PIVOT i ICOLUM)* IPIVOTI ICOLUM >f l 

INTERCHANGE ROWS TO PUT PIVOT ELEMENT 
IF CIROW-ICOLUM) 140* 260* 140 

140 deter**-deterh 

DO 200 L* 1 *N 
SWAP* A ( I ROW « L ) 

A<1 WOW • L I * A I lCOLUM*L ) 

200 A < ICOLUM*L I *S*<AP 

IF(H> 260 • 260* 210 
210 00 250 L*I * M 
SWAP«B( IRON *L ) 

B( I WOW * L ) *B ( ICOLUM.L ) 

250 QUCOUUM.l )*SWAP 
260 1 wK <1*1 ) = I ROW 

1 WKC 1 *2 > - I COLUM 
P 1 VOT » A < I COLUM * 1 COLUM 1 
IF ( rOP*EQ*l ) GO TO 321 
IF (PIVOT 1 I00u*l 06* I OOC 

SCALE the DETERMINANT 

c 

luou P I VOT I = P 1 VOT 

If <Ae*S<DETE«M)-Rl >10 30*1010*1010 
1 01 l DE rtiRM = DE fERM/Rl 

: scale * I scale* i 

IF CABS(QLT£«M)-Wi >1060# 1020* 1020 
1^2 0 OETi:RM*DETERM/Rl 
I SCALE* ISCALE-M 
GG TO 1 060 

1 ^3t> I F ( AQS ( OETtRM )-R2 >10*0*1040*1 060 
1C4Q OKTERMsDETERMARi 

ISCALC=1SCALE-1 

IF (At3S(OETERM)-«2 >1050* 105 l* 1060 

1050 determ=detewm*ri 
I SCALE* I5CALE-1 

1 v>60 IF (AfiS( PI VOT l ) -W 1 ) 1 090* 107 0* 1 0 70 
1070 P!V0TI=PIV0TI/R1 
I SCALE* ISCAlE-M 

IF (AdSIPl VOT I )-R| >320* 1 080* 1 OBO 
1080 PI VOT I *P1 VOT I/Rl 
I SCALE* 1 SCAlE+ 1 
GO TO 320 

1 090 I F ( AOS \ P 1 VOT 1 )-«2 ) 2000 * 2000 * 320 


ORIGINAL PAG 3 IS 
DP POOR QUALITY 


ON 01 AGONAL 


49 



APPENDIX D 


2OU0 PIVOT! »PI VOT l *Wl 
1 SCALE * 1 SCALE* i 

|F(A8S(PIV0TI >-R2 >2010*2010*320 
2wl0 P ! VOT I «P l VOT I *P l 
l SCALE »i SCALE- l 

320 DETt«M«OCTERM*PlvOTI 
C 

C 01VIDE PIVOT ROW BY PIVOT ELEMENT 

C 

32 1 IF <P! VOT >330* 106.330 
330 A ( ICOLUM* ICOLUM >* i .0 

DO 350 L l 1 *N 

350 A< ICOEUM.L >»A< ICOLUM.L >/PI VOT 
IFIM) 380* 360 « 360 
36u DO 370 L= I *M 

37. B( t COLUM *L ) *B < I COLUM *L ) /P l VOT 
C 

C WFDUCE NON-PIVOT WO VS 

C 

380 DO 550 L I * 1 *N 

IF (L l - ! COLUM > 400* 550* 400 
400 T»A(LI . ICOLUM) 

A ( L l . ICOLUM)«0*0 
DO 450 LM • N 

450 A(L1«L)*A<LI*L) -A < 1 COLUM *L ) * T 
IF(M) 550. 550. 460 
460 DO 500 L s 1 *M 

600 b(Ll *L>*BtLl *L >-B ( ICOLUM.L >*T 
550 CONTINUE 
C. 

C INTERCHANGE columns 

c 

DO 710 1*1 *N 
L«N+ l - I 

IF < IWKIL.i )-l*K<L«2 > >630*710*630 
630 JWOW* 1 WK ( L * 1 > 

JCOLUM* I VK ( L • E ) 

DO 705 K*1 * N 
SWAP* A (K * JWOW > 

A <K* JWQV ) * A (N* JCOLUM ) 

A (K « JCOLUM )*SWAP 
7C5 CONTINUE 
710 CONTINUE 
**40 RETURN 
ENO 
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TABLE I.- SIMtoHY OF BOUNDARY CONDITION EQUATIONS 3 


Collective node 

Cyclic node 

Scissors mode (no 
hinge offsets) 

Scissors mode 
(hinge offsets) 

L - <K op • Hop'*) 2 ) 6y - 0 

6y - 0 

6y - 0 

b (6y * 0) 

6 x “ 0 

D - (K ip - M ip 4> 2 ) 6 X ■ 0 

6 X - ° 

C (6 X - 0 ) 

8*0 

M - Kg8 - 0 

8 * C 

b (M - Kg 6 * cn 

Q - Rpt * 0 

i|> * o 

4* * o 

c (C - * 0) 

d (T - K<4» * 0) 

d {T - Kc<J> « 0) 

d (T - K c 4> - 0) 
®<4» * 0) 

d (T - K c 4> = 0) 


a Boundary condition applied at center of rotation (station 1) unless other- 
wise noted. 

^Boundary condition applied a: flapping hinge offset. 
c Boundary condition applied at lagging hinge offset. 

^Boundary condition applied at pitch horn offset. 

e For two-bladed rotors only and applied at :.itch horn offset. 
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TABLE II.- DESCRIPTION OT PROGRAM INPUT REQUIREMENTS 


Card 

POtIMt 

Variabla 

Units 

Dascr ip t ion 

typ*| 





i 

«(A4,6X) 



Program opt.onsi DECK. NAME. PUN?, 





SNAP, MODE. ALLM. PLOT. TORS, and 





ULSW 

2 

It,4X, 10A« 

NR Ml' 


Problan nuabar 



ITLE 


Problan tit la 

3 

•PlO.O 

XT (X T ) 

in-lb/rad 

Mast torsional ctiftnass par blada 



MOP (Mop) 

lb-sacVin. 

NUb oot-of-plana inart ia par blada 



M. (M|p) 

lt-sac 2 /in. 

Nub in-plana inart ia par blada 



KOP (Hop) 

lb/in. 

Rub out -of -pi ana stiffnass par blada 



XIP (K^p) 

lb/in. 

Hub in-plan* stiffness per blade 



XBBTA (Kg) 

in-lb/rad 

Blada flapping spring 



KPSI (XT) 

in-lb/rad 

Blada lagging spring 



KC (K c ) 

in-lb/rsd 

Control systam stiffnass 

4 

2lS # tPS,0,F10.0 

N 


N unbar of sag nan ts 



JHUB 


Nuabet of nonfaathar ing sagnants 

t 



AZBAk (I) 

in. 

f 

Segment length for uniform 




, 

distribution ! 


1 

TWIST <0 t ) 

dag 

Linear twist from center of fetation j 


i 

1 



to tip | 


| 

BODES 


i 

Number of blades j 



cuorr 

in. 

Lading binge offset from center lit.*?! 

1 



rHOFT 

ir*. 

flapping hinge offset from center * 





line 


1 

i 

i 

PHOFF 

in. 

Titch horn off* w fns centei lint 

WV, 

1 
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TABLE It.- Concluded 


Card 

type 

Format 

Variable 

Unite 

Description 



BLAST 

c*» 

Upper limit of frequency calculations 

5 

I5.10F5.0 

IBM 


Number of rotor speeds 



MUM 

tpm 

Rotor speeds 

• 

15, 10P5.0 

1HCOL 


Number of rotor collective pitch 
angles 



BCQLL <e c ) 

de9 

Rotor collective pitch angles 

? 

8F10.0 

i 

•l 

in. 

Outboard end stations of segments 

9 

8F10.0 

a KIB (EI b ) 

lb-in 2 

Beamwise bending stiffness 

9 

8F10.0 

a RIC (EI C ) 

lb-in 2 

Chordwise bending stiffness 

10 

8F10.0 

tVPPL (w,N) 

lb/in, lb 

Weight per unit length and tip weight 

11 

8F10.0 

C THD (6) 

deg 

Nonlinear twist distribution at 
station 

12 

8F10.0 

1 

•BYEB (1^) 

lb- sec 2 

Beanwise second mass moment of 

inertia j 

| 13 

1 

i 

8F10.0 

| 

*EYEC (1^.) 

lb-aec 2 

Chcrdwise second mass moment of i 

inertia ! 

i 

14 

j 8F10.0 

*OJ 

lb- in 2 

! 

Torsional stiffness , 

15 

8F10.0 

a SB (S b ) 

in. 

Beamwise shear center offset 

i 

16 

8F10.0 

a SC <S C ) 

in. 

Chordwise shear center offset 

i 

17 

8F10.0 

d RB (r b ) 

in. 

Beamwise offset of center of gravity 

18 

8F10.0 

d RC (r c ) 

in. 

Chordwise offset of center of gtavity 


a N values are input, one for each segment. 

&N values are input and the tip weight is listed in field N ■* 1 . 

C N ♦ 1 values are input, one for each station. 

values are input and the offsets pertaining to the tip weight are 
input in field N ♦ 1 . 
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cud tn* 


10 


12 


13 


14 


15 


U 


17 


18 


C*rd col tan j 

looopooooon I 111 1 1 1 1 222222222233333333334444444* 

4*6555655555666666666677777777778 | 

12345678901234567890123456789012345678901234567890123456789012345678901234567890 

DFCK 

TORSION PLOT *Of> SMAPF 


PI INCH 



100000 

COMPtHFR PROGRAM SA*PlF C*SF 



I 

166000* 

.99 .269 1160. 1150. 


0. 

0. 

16500. 

30 2 

.0 .0 4. 3. 3. 3. 16000. 





3 1100 

1200 1300 





3 0 

6 12 





1.5 

3.0 4.935 6.87 9 

.1 

11.235 

13.37 

14.87 

17.1 

18.8 20.5 22.5 25. 

179 

27.857 

30.536 

33.214 

35.893 

38.571 41.26 43.929 46. 

607 

49.286 

51.964 

54.643 

ST. 321 

60. 61. 61.25 61. 5 

62.05 



100000000. 

100000000.1500000. 1500000. 500000. 


182900. 

182900, 

1198000. 

500000. 

135000. 135000. 79400. 62900. 


62900. 

62900. 

62900. 

62900. 

62900. 62900. 62900. 62900. 


62900. 

62900. 

62900. 

62900. 

62900. 100000. 100000. 10000. 


10000. 



100000000. 

100000000.1500000. 1500000. 1250000 


911400. 

911400. 

11 98000. 

1250000. 

972000. 972000. 820000. 560000, 


560000. 

560000. 

560000. 

560000. 

560000. 560000. 560000. 560000. 


560000. 

560000. 

660000. 

560000. 

560000. 1000000. 1000000. 100000. 


1 oOOOO. 



1.7083 

1.7083 .4185 .4185 .1804 

.0705 

.0705 

. 1226R 

.1652 

.02777 .02777 .02779 .02766 

.02766 

.02766 

.02766 

.02766 

.02766 .02766 .02766 .02766 

.02766 

.02766 

.02766 

.02766 

.02766 .04315 .08273 .03744 

.01101 

0.0 


0. 

0. 0. 0. 

0. 

0. 

0. 

0. 

0. 

0. 0. 0. 

0. 

0. 

0. 

0. 

0. 

0. 0. 0. 

0. 

0. 

0. 

0. 

0. 

0. 0. 0. 

0. 

0. 



.001 

.001 .00057 .00057 .000143 

.0000285 

.0000286 

.000062 

.000114 

.000023 .000023 .000064 .000063 

.0000*3 

.00006* 

.000063 

.000063 

.000063 .000063 .000063 .000063 

.000063 

.000063 

.000063 

.000063 

.000063 .000131 .000074 .000068 

.000023 



• 1 F ♦ 1 0 

.ie+101000000. 1000000. 2!3SOO. 


202000. 

202000. 

910500. 

213500. 

109400. 91390. 75900. 75900. 


75900. 

76900. 

75900. 

75900. 

75900. 75900. 75900. 75900. 


75900. 

75900. 

75900. 

75900. 

75900. 114000. 114000. 11400. 


11400. 



0. 

0. 0. 0. 

0. 

0. 

0. 

0. 

0. 

0 . 0 . 0 . 

0. 

0. 

0. 

0. 

0. 

0. 0. 0. 

0. 

0. 

0. 

0. 

0. 

0. 0. 0. 

0. 

0. 



0. 

0. 0. 0. 

0. 

0. 

0. 

0. 

0. 

0 ♦ 0 . 0 . 

0. 

0. 

0. 

0. 

0. 

C. 0. 0. 

0. 

0. 

0. 

0. 

0. 

0. 0. 0. 

0. 

0 . 



0. 

0 . 0 . 0 . 

0. 

0. 

0. 

0. 

0. 

0. 0. 0. 

0. 

0. 

0. 

0. 

0. 

0. 0. 0. 

0. 

0. 

0. 

0. 

0. 

0. 0. 0. 

0. 

0. 

0.0 


0. 
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TABLE VII.- DESCRIPTION OP PUNCHED OCTPUT DATA 


Card 

Description 

F*->n«at 

1 

Identification card containing case nuaber and 
title 

16, 4X, 10 A6 

2-4 

Distribution of weight per unit length and tip 
weight, lb/in. and lb 

7F10.5 

5-7 

Beaswise second Mass aoawnt of inertia, lb-sec? 

7P10.5 

8-10 

Chordwise second Mass ncasent of inertia, lb-sec? 

7F10.5 

11-76 

Collective node deflection shapes, case number, 
■ode ninber. Mode type 

6F10.6, 10*, 16, 212 

77-142 

Cyclic mode deflection shapes, case number # mode 
number, and mod* type 

6F10.6, 10*, 16, 212 

143-208 

Scissors mode deflection shapes, case number, 
mode number, and mode type 

6F10.6, 10X, 16, 212 

209-214 

Collective mode cyclic detuni, ; information, ref- 
erence values of collective Pav^*. and rotor 
speed, case number, mode number, and mode type 

6F10.2, 2F5.0, 16, 212 

215-220 

| 

Cyclic mode cyclic detuning information, refer- 
ence values of collective pitch and rotor 
speed, case number, mode number, and mode type 

6FU.2, 2F5.0, 16, 212 

221-226 

Scissors mode cyclic detuning information, refer- 
ence values of collective pitch and rotor 
speed, case number, mode number, and mode type 

6F10.2, 2F5.0, 16, 212 
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Original page is 

oe POOR QUALITY 


nan vm.~ ttumuarm or nogdm optional output ( 


CA» MCI) 



Card coluan 

Card no. 

*000000401 11111111 l??????????313333?7?:H4*444*44***5*« , 5**5*A**66f 444' *7777777777* 


1 ??**67O«01 33***7*901 73***70901 ?3*©f7*90 173***7990 I3?**67*90l?3**f 7*001*34**7*0* 

1 

190000 

rr*puTF9 ponrio** c**f 





2 

1.ft***9 

.41**0 

.77*73 

.074*0 

.11*13* 

.09*49 

. 0??7A 



3 

,037*0 

.037** 

.07746 

,07744 

.037©* 

.0776© 

.037*6 



4 

.037** 


.07744 

,07744 

.4 7746 

.07493 

.40064 



5 

o.r,onpo 

6, 66464 

n.40000 

0.00400 

0 . 004*0 

4.00000 

0.04400 



C 

o.ocooo 

O.ooncr 

0.00000 

0.00400 

n.nnono 

0. 00044 

0.44600 



7 

0.00000 

o.*oooo 

0*00000 

0 .pooflft 

0.00444 

0.00060 

o.oooon 



• 

.000*9 


.000?3 

.00403 

.04004 

.00007 

.0404* 



» 

.00004 

. 00006 

.00004 

.00004 

.00046 

,00006 

.04446 



10 

,00*04 

.00004 

.90006 

,00004 

.400*4 

.0*00* 

.04000 



143 

0.004060 

n.ononno 

-.oonooo 

-.004400 

.noni*? 

-.004000 

1 onnoo 

1 

0 

144 

•.000004 

.00**37 

-.rooooo 

-.000016 

.40*737 

-.004 000 

1 00040 

1 

4 

145 

-.000071 

,r.I30*» 

-.000000 

-.00i'0?4 

.4173** 

-.000004 

1 40660 

1 

6 

144 

-.00007© 

.07140* 

-.ooooon 

-.ocnnsi 

.4760*0 

-.00004* 

164600 

1 

0 

147 

-,0000'M 

.03030* 

-.000000 

-.00447* 

.93*761 

-.000000 

104660 

1 

0 

144 

-.00007* 

.039143 

-.000000 

-.000471 

.4*3*36 

-.000000 

1 00664 

1 

0 

149 

-.000014 

.P*T93*> 

-.oooooo 

-.444011 

.4*73*7 

-.000000 

100644 

1 

n 

150 

-.000004 

,49*7** 

-.000000 

.000400 

.0611*6 

-.000000 

1 44466 

1 

4 

151 

.000004 

,99 C M? 

-.000000 

.44041? 

.070040 

-.ooooon 

1 60440 

1 

0 

152 

,00001© 

,074*71 

-.000000 

.00<19?e 

.07*90? 

-.000000 

104400 

1 

4 

153 

.000071 

.P9331-* 

-.000000 

.7947*7 

0,400400 

.920009 

1 60000 

1 

0 

154 

0.000000 

0.000000 

-.000091 

.000139 

-.600000 

-.000001 

104000 

7 

4 

155 

.004340 

-.00000* 

-.000001 

.00**47 

-.000016 

-.000601 

106006 

? 

6 

154 

.oirm 

-.000077 

-.400001 

.01703? 

-.400077 

-.600001 

1 60646 

7 

6 

157 

.nann? 

-.000031 

-.000001 

.07**0? 

-.004037 

-.000001 

106660 

? 

4 

154 

.079937 

-.000037 

-.000, '01 

.03*303 

-.040079 

-.000001 

160460 

7 

6 

159 

.03*497 

-.00007* 

-.000001 

• 0*3110 

-.44001 P 

-.064001 

1 00646 

7 

4 

140 

.0*793© 

-.00001? 

-.000001 

.0519*7 

-.onooo© 

-.000001 

100000 

7 

0 

141 

.0*4437 

.000003 

-.000001 

.060903 

.000017 

-.000001 

100006 

7 

6 

142 

.04*370 

.000071 

-.000001 

,069*4? 

.000031 

-.000041 

146060 

? 

0 

143 

.0743*0 

.000041 

-.000001 

.07***! 

.0000*1 

-.ormoi 

144064 

? 

6 

144 

.043733 

.00004? 

-.000001 

1 ,0*?OQ© 

0 , ooooon 

.470044 

1 00600 

? 

6 

145 

o.oooooo 

0.000000 

-.ooooon 

-.0on?o© 

-.0*0037 

-.400000 

1 60066 

3 

4 

144 

-.0093?* 

-.0011*4 

-.0*0000 

-.01*1©1 

-.0077*1 

-.oooorn 

1 6446r 

7 

6 

147 

-.07470* 

-.003740 

-.00*000 

-.033740 

-.00*107 

-.004440 

140046 

7 

n 

144 

-.039*90 

-.00**77 

-.000000 

-.0*41*0 

-.00**70 

-.000064 

166004 

3 

4 

149 

-.0*4*1* 

-.00*70© 

-.000000 

-.0*47*3 

-.00*6*6 

-.000044 

1 66466 

3 

6 

170 

-.0*3*4* 

-.00*3*9 

-.000000 

-.03*67? 

-.00*7*0 

-.000040 

1 00640 

7 

6 

171 

-.0314*1 

-.0*3*70 

-.oooroo 

-.07763* 

-.0077*7 

-.000000 

106466 

3 

6 

12 

-.01147! 

-.non** 

-.OOCOOO 

.001714 

.0401 ©6 

-.004004 

tooooo 

7 

6 

173 

.01**03 

.0019*4 

-.rooooo 

.0314*4 

.007974 

-.004044 

1 66606 

7 

4 

174 

,0***7* 

.00* c 70 

-.ponton 

,0^o 04 

.00*077 

-.000404 

1 Conor 


6 

175 

.0*333? 

.010707 

-.000000 

7. *70757 

0.004040 

, o?44no 

J 64066 

7 

4 

174 

0.000000 

o.oonnno 

-.ooonoo 

,0000*0 

-.00033? 

-.400046 

1 06666 

k 

6 

177 

.0017*7 

-.010739 

-.000000 

.007*0© 

-.0197*3 

-.444400 

144466 

k 

6 
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MU VI XI.- Concluded 


Cirt no. 

Ctrl eoium 

IN 

.003*91 

•.024273 

-.000000 

,004347 

-.034456 

-.oooooo 

1 OftftftO 

4 

0 

m 

.0090.0 

-.041507 

-.000000 

.00**39 

-.04*415 

— .oooooo 

100000 

4 

0 

in 


-.044201 

-.000000 

,00**46 

-.04741? 

-, OOOOOO 

tonoeo 

4 

0 

in 

,«MWt 

•.0**773 

—.000004 

,004142 

-.03*179 

-.oooooo 

100000 

4 

0 

112 

.143*93 

*-.031999 

-.090000 

,002420 

-.021327 

-.oooooo 

100OO0 

4 

0 

113 

.201235 

•.004*24 

—.000000 

-.000411 

,003*34 

-.oooooo 

lOOOOO 

6 

0 

1M 

-.0*2234 

.01*03* 

-.00O000 

-.004132 

.034127 

-.oooooo 

1 00«0 o 

4 

ft 

IIS 

-.00*901 

.0*0210 

— . 000000 

-,00*?22 

.022246 

-.oooooo 

1 onooo 

4 

0 

IH 

••010399 

,00333“* 

oooooo 

3. *9907* 

n.onoooo 

,020000 

looono 

4 

0 

117 

0.000000 

0.000000 

—.000001 

,000397 

.000032 

-.000001 

100000 

4 

0 

IN 

.011004 

.001141 

-.000001 

.0910*9 

.0099*5 

-•0000O1 

looooo 

4 

0 

Ilf 

.027524 

.00993* 

-.oooooi 

,030200 

.00391* 

-.oooooi 

100000 

4 

ft 

IN 

.030277 

.003*1* 

-.000001 

,09*7*1 

.00261* 

-.000001 

1 00000 

e 

0 

111 

.013030 

.0OU07 

-,040001 

-,002224 

-.00026* 

00*001 

1 OOftftft 

c 

ft 

112 

•*01171* 

-.002124 

— , oooooi 

-.032044 

-.003*67 

-.oooooi 

100000 

5 

ft 

1t3 

-.041007 

-.0099** 

-, 000001 

-.05**01 

-.00604* 

-.000001 

100000 

4 

0 

m 

-.097310 

-.000115 

-, 009001 

-.040164 

-.004346 

-.oooooi 

1 OftftOO 

4 

0 

Iff 

034221 

-.003030 

-.000001 

-.010420 

-.001144 

-.OOOOOI 

1 0000ft 

4 

ft 

m 

.017701 

.*01414 

-.OOOOOI 

.0**001 

,005344 

-.oooooi 

lonooft 

* 

ft 

117 

.003333 

.00441* 

-.000001 

4.214033 

0,000000 

,020000 

1 Oftftftft 

« 

ft 

IN 

o.nooooo 

0.000000 

10.9*U«6 

•002421 

-.010704 

10,222430 

100000 

6 

0 

IN 

.©0?7*4 

-.012437 

10,700*79 

.002272 

-.0141*2 

12.2410** 

100000 

2 

ft 

2N 

.00*33? 

-.013440 

1*. 1*773* 

.00120* 

-.01 0**2 

15.29*204 

IftOftOft 

6 

ft 

Ml 

•.000213 

-.0047*0 

17,019179 

-.001262 

-.001214 

22.042272 

inoooo 

6 

0 

M2 

•*003440 

.00*244 

90,09099* 

-.004554 

,0120*6 

31.055474 

100OO0 

6 

0 

M3 

•.004100 

.014255 

39,19*910 

-.004423 

.023033 

39.063002 

looooo 

6 

ft 

204 

-.003073 

.025754 

*9.09710* 

-.002076 

,02593? 

45.261*41 

1OOO00 

6 

ft 

MS 

•.000440 

.023447 

*0.7*30*0 

.006477 

•01*174 

*1.24415? 

1 00000 

2 

0 

2N 

.001041 

.01014* 

«“». 31*931 

,002036 

• 0001 13 

44,464**6 

Iftftoro 

6 

0 

M7 

.001010 

-.011347 

**.19*7*9 

,000276 

-.023444 

*2 ,4**32 1 

inoooo 

* 

ft 

Ml 

..000241 

-.034234 

«*. 9*9*00 

7.460446 

o, nooooo 

,020000 

10000ft 

2 

0 

221 

320.40 

324.44 

309,70 

325.47 

6,00 

100,00 

6.1200,100000 

1 

0 

222 

1140.34 

1144.24 

139*. 79 

13*4,62 

6.00 

100.00 

6.1200.10000ft 

? 

0 

223 

2027.10 

2414.20 

399?. *9 

3234,12 

6,00 

100,00 

2.1200, lOOOftft 

3 

0 

224 

40*0.74 

4044,72 

*300.07 

4317.64 

6,00 

100,00 

2,1 200 , 1 OOOftft 

4 

ft 

223 

**20.00 

**21.44 

*030,11 

6032.1? 

6.00 

1 on, on 

^.1 200, 1 OOOftft 

« 

0 

*; s 

4*42,00 

«*?7,to 

»**T,7* 

4460,10 

6,00 

100,00 

6.1200,100000 

6 

ft 

_ J 
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TABLE IX.- UNIFORM BEAM ANALYTICAL REPRESENTATIONS USED 


FOR CORRELATION STUDIES 


Property 

Case 1 

Case 2 

Case 3*~ 

Case 4* 

<K) b , lb-in 2 .... 

8000 

8000 

8000 

8000 

(£l) c , lb-in 2 .... 

32 * 10* 

32 x 10 9 

32 X 10 9 

32 x 10 9 

GJ, lb-in 2 

8717 

8717 

8717 

8717 

w, lb/in. 

0.024 

0.024 

0.024 

0.024 

W, lb 

0 

0*6 

0 

0 

!&>* lb-sec 2 

-0.5176 x 10-5 

-0.5176 x 1 0“5 

-0.5176 x 10-5 

-0.5176 x 10-5 

I Ac, lb-sec 2 

0.1087 x 10- 3 

0.1087 x 10*3 

0.1087 x 10-3 

0.1087 x 10-3 

N 

25 

25 

25 

25 

i, in 

1.0 

1.0 

1.0 

1.0 

S b , in 

0 

0 

0 

0 

S c# in* ..«•••• 

0 

0 

0 

-0*7454 

r b . in 

0 

0 

0 

0 

rpt in. ....... 

0 

0 

0 

0 

ft, rpm 

0 

0 

0 

0 

0 C » 

0 

0 

0 

0 

*»9 

0 

0 

0 

0 

Kc, in-lb/rad .... 

0 

0 

CL 

OP 

Rip, lb/in 

0 

0 

<x> 

ft' 

Rop, lb/ in. ..... 

0 

0 

0 

0 

Rj» in-lb/rad .... 

0 

0 

ft* 

ft' 

Rg, in-lb/rad .... 

0 

0 

0 

G 

R^, in-lb/rad .... 

0 

0 

0 

0 

Mip, lb- sec 2 / in. . . . 

0 

0 

0 

0 

Mop, lb-sec 2 / in. . . . 

0 

0 

0 

. . J 

0 


•Program was modified to achieve proper boundary conditions at the outboard 

tip. 
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TABLE X.- ROTATING PROPELLER ANALYTICAL REPRESENTATION USED 


FOR CORRELATION STUDIES (CASE 5) 


Z, 

in. 

(EI) b , 

lb- in 2 



(EI)c, 
lb- in 2 

GJ + EBi 0' 2 , 
lb-in 2 

w# 

lb/ in. 

deg 

J bb« - 
lb-sec^ 

*OC'« 

lb-sec* 

4.0 

OO 

OO 

CO 


0 



5.0 

10 000 000 

300 000 000 

1 000 000 

0.5000 

-2.55 

0.000010000 

0.002000 

6.0 

800 000 

100 000 000 

700 000 

.3993 

-5.10 

.000007137 

.001692 

7.0 

180 000 

51 000 000 

575 000 

.3100 

-7.65 

.000005069 

.001400 

8.0 

140 000 

47 900 000 

480 000 

.2761 

-10.20 

.000003001 

.001265 

9.0 

110 000 

46 800 000 

420 000 

2600 

-12.75 

.000002507 

.001225 

10.0 

92 000 

45 800 000 

382 000 

.2437 

-15.50 

.000002013 

.001184 

n.o 

72 000 

45 100 000 

348 000 

.2408 

-18.25 

.000001755 

.001204 

12.0 

63 000 

44 400 000 

320 000 

.2379 

-21.00 

.000001496 

.001223 

13.0 

53 000 

44 100 000 

295 000 

.2294 

-23.55 

.000001270 

.001228 

14.0 

42 000 

43 800 000 

275 000 

.2209 

-26.05 

.000001043 

.001233 

15.0 

37 000 

43 800 000 

261 000 

.2166 

-28.45 

.000001042 

.001245 

16.0 

32 000 

43 800 000 

251 000 

.2123 

-31.05 

.000001042 

.001256 

17.0 

29 500 

44 100 000 

241 000 

.2081 

-33.45 

.000000882 

.001263 

18.0 

27 000 

44 400 000 

230 000 

.2039 

-35.55 

.000000722 

.001270 

19.0 

26 500 

45 800 000 

214 000 

.2018 

-37.60 

.000000691 

.001289 

20.0 

26 000 

47 200 000 

196 000 

.1997 

-39.55 

. 00000065° 

.001307 

21.0 

25 500 

49 300 000 

173 000 

.1955 

-41.50 

.000000638 

.001311 

22.0 

25 000 

51 400 000 

153 000 

.1912 

-43.00 

.000000616 

.001314 

23.0 

24 500 

53 700 000 

125 000 

.1912 

-44.30 

.000000627 

.001340 

23.5 

24 000 

54 900 000 

100 000 

.1912 

-44.90 

.000000632 

.001353 

24.0 

24 000 

56 000 000 

u 

78 000 

.1912 

-45.55 

.000000638 

.001366 
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TABU XI.- COMPARISONS OF COMPUTED FREE-FREE BENDING NATURAL FREQUENCIES AND 
MODE SHAPES WITH EXACT SOLUTIONS FOR A NONROTATING, UNIFORM BEAM (CASE 1) 


Natural frequency, w, cpm 


First soda 
[Exact 


| Computed 
ME. 68 


ME. 88 


Second aoda 

Third aoda 

Fourth 

aoda 

Fifth mode ] 

Computed 

Exact 

Computed 

Exact 

Computed 

Exact 

Computed 

1 

Exact ■ 

2667.9 

2673.5 

5225.8 

5241.2 

8631.3 

8663.9 

12883 

12942 


Node shape 


L> 

First mode 

Saeond mode 

Third inode 

Fourth n 

tode 

Fifth mode 

Computed 

Exact 

Computed 

Exact 

Computed 

Exact 

Computed 

Exact 

Computed 

**xact 

0 

-0*608 

-0. 608 

o 1 

H 

0 

0.714 

0.711 

0 

0 

-0.714 

-0.707 

.04 

-.605 

-a 604 

-.108 

-.108 

a 697 

.694 

. 199 

.198 

-.672 

-.665 

.08 

-.594 

-.594 

-.212 

-.212 

.647 

• 6%4 

.381 

♦ 379 

-.551 

-.545 

.12 

-.576 

-.576 

-.312 

-.311 

. 566 

.564 

.534 

.531 

-.365 

-.360 

.16 

-.552 

-.551 

-.403 

-.403 

a 459 

.457 

.645 

.641 

-.136 

-.132 

.20 

-.521 

-.520 

-.484 

-.483 

.330 

. 128 

.705 

.700 

.109 

.111 , 

.24 

-.483 

-.482 

-.551 

-.551 

.186 

.184 

.709 

.704 

.342 

.341 i 

.26 

-.439 

-.438 

-.605 

-.604 

.034 

.032 

.657 

.652 

.534 

.531 | 

.32 

-.389 

-.389 

-.642 

-.641 

-.119 

-.120 

* 555 

.549 

.664 

.659| 

.36 

-.333 

-.333 

-.662 

-.661 

-.265 

-.266 

.408 

.403 

.716 

.709] 

.40 

-.272 

-.272 

-.663 

-.662 

-.397 

-.397 

.231 

.226 

• 684 

♦ 675 | 

.44 

-.207 

-.206 

-.646 

-645 

-.508 

- 507 

.036 

.032 

.572 

.563 

.48 

-.136 

-.136 

-.611 

-.609 

-.592 

-.591 

-.160 

-. 163 

.394 

.384 

.52 

-.062 

-.062 

-. S57 

-.556 

-.645 

-.643 

-.341 

-.343 i 

.171 

.162 

.56 

.016 

.016 

-.486 

-.485 

-.663 

-.661 

-.493 

-.493 

-.070 

-.077 

.60 

.097 

.098 

-.399 

-.397 

-.646 ! 

-.643 | 

-.602 

-.600 1 

-.301 

-. 305 i 

.64 

. 161 

.182 

-.297 

-.295 

-.592 

-.588 

-.658 

-.655 

] -.492 

| -.493 i 

.68 

.268 

.268! 

-.182 

-.180 

-.503 

-.499 

-.657 

-.652 

I -.620 

| -.617| 

.72 

.356 

.356 

-.055 

-.053 

-.381 

-.377 

-.595 

-.588 

i -.667 

i -.661 , 

.76 

.446 1 

.446 

.082 

.084 

-.231 

-.226 

-.474 

-.467 

-.625 

1 -.616 ; 

.80 

.537 

.537 

.226 

.227 

-.056 

j-. 052 ! 

-.302 

-.294 

-.494 

! -.483 1 

.84 

.629 1 

.629 

.375 

.377 

.137 

! .142 

-.086 

-.078 

' -. 284 

1 -.272 ; 

.88 

.721 

.721 

.52« 

.530 

.344 

.347 

J .163 

i .170 

-.010 

.001 

.92 

.814 

.814 

.685 

.686 

.559 

.562 

.433 

! .438 

.308 

.317 

.96 

.907 

.907] 

.842 ! 

.843 

.779 I 

.786 

.715 

.718 

.650 

.655 

1.00 

1.000 

1.000] 

1.000 

i 

1.000 

1 

1.000 

1.000 

1 

1.000 

! 1.000 

1 j 

1.000 

1. 000 





TABLE XII.- COMPARISONS OF COMPUTED CLAMPED- FUSE BENDING NATURAL FREQUENCIES AND 
NODE SHAPES NITH EXACT SOLUTIONS FOR A NONNOTATING, UNIFORM BEAM (CASE 1 ) 































TABLE XIII.- COMPARISONS OP COMPUTED FREE-FREE TORSION NATURAL FREQUENCIES AND 
MODE SHAPES NITP EXACT SOLUTIONS FOR A NONROTATING, UNIFORM BEAM (CASE 1) 


Natural frequency, <*■’ , cpm 


First node 

Second node 

Third node 

Fourth node 

Fifth node 

Computed 

Exact 

Computed 

Exact 

Computed 

Bx»ct 

Computed 




11004 

11012 

21965 

22023 

32840 

33035 



54157 

55058 


Mode shape 


«j/R 

iirst node 

Second k,Jde 

Third node 

Fourth node 

— 

Fifth node 

Computed 

Exact 

Computed 

m 


Exact 

Computed 



Exact 

0 

-1.000 

-1.000 

1.000 

1.000 

- 1.000 

-1.000 

1.000 

1.000 

- 1.000 

-1.000 

.04 

-.992 

-.992 

.969 

.969 

-.930 

-.930 

.876 

.876 

-.809 

-.809 

.08 

-.969 

-.969 

.876 

.876 

-.729 

-.729 

.536 

.536 

-.309 

-.309 

.12 

-.930 

-.930 

.729 

.729 

-.426 

-.426 

.063 

.063 

.309 

.309 

. 16 

-.876 

-.876 

.536 

.536 

-.063 

-.063 

-.426 

-.426 

.809 

.809 

.20 

-.809 

-.809 

.309 

.309 

.309 

.309 

-.809 

• 809 

1.000 

1.000 

.24 

-.729 

-.729 

.063 

.063 

.*17 

.637 

-.992 

-.992 

.809 

.809 

.28 

-.637 

-.637 

-.187 

-.187 

.876 

.876 

-.930 

-. 93C 

.309 

.309 

.32 

-.536 

-.536 

-.426 

-.426 

.992 

.992 

-.637 

-.637 

-.309 

-.309 

.36 

-.426 

-.426 

-.637 

-.637 

.969 

.969 

-.187 

-.187 

-.809 

-.809 

.40 

-.309 

-.309 

-.809 

-.809 

.809 

.809 

.309 

.309 

- 1.000 

- 1.000 

.44 

-.187 

-.187 

-.930 

-.930 

.536 

.536 

.729 

.729 

-.809 

-.809 

.48 

-. Of 3 

-.063 

-.992 

-.992 

.187 

.187 

.969 

.969 

-.309 

-.309 

.52 

.063 

.063 

-.992 

-.992 

-.187 

187 

.969 

.969 

.309 

.309 

.56 

187 

.187 

-.930 

-.930 

-.536 

-.536 

.729 

.729 

.809 

.809 

.60 

.309 

.309 

-.609 

-.809 

-.809 

-.809 

.309 

.309 

1.000 

1.000 

.64 

.426 

.426 

-.637 

-.637 

-.969 

-.969 

| -.187 

-. 187 

.609 

! .809! 

.68 

.536 

.536 

-.426 

-.426 

-.9.2 

-.992 

1 -.637 

-.63/ 

.309 

i . 309 : 

.72 

.637 

.637 

-.187 

-.187 

-.876 

-.876 

, -.930 

-.930 

-.309 

| -. 309! 

.76 

.729 

.729 

.063 

.063 

-.637 

-.637 

-.992 

-.992 

-.809 

-.809 

.80 

.809 

.809 

.309 

.309 

-.309 

1 -.309 

! -.809 

-.809 

- 1.000 

- 1.000 

.84 

.876 

.876 

.536 

.536 

.063 

; .063 

-.426 

-.426 

-.809 

-,809 

.88 

.930 

.930 

.729 

.729 

.42b 

.426 

.063 

.063 

-. 309 

j -.309 

.92 

.969 

.969 

.876 

.876 

.729 

.729 

.536 

j .536 

.309 

j .309 

.96 

.992 

.992 

.969 

.969 

.930 

.930 

.876 

.876 

.609 

.809 j 

1.00 

1.000 

1.000 

1.000 

’.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 
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TABLE XIV.- COMPARISONS OP COMPUTED CLAMPED- FREE TORSION NATURAL FREQUENCIES AM) 
NODE SHAPES WITH E3UCT SOLUTIONS FOR A NONROTATING, UNIFORM BEAM (CASE 1) 

Natural frequency, u, epa 


First aode 

Second bode 

Third node 

Fourth aode 

Fifth aode 

C aputed 

Exact 

Computed 

Exact 

Computed 

Exact 

Coaputed 

Exact 

Coaputed 

Exact 

S5C4.9 

S505.8 

16493 

16517 

27416 

27529 

38231 

38541 

48895 

49552 


Node shape 


Xj/R 

First 

aode 

Second 

aode 

Third node 

Fourth aode 

Fifth node 

Coaputed 

Exact 


Exact 

Coaputed 

Exact 

Coaputed 

r 

Exact 

Coaputed 

Exact 

0 

0 

0 

0 

0 

mam 

0 

0 

0 

0 

0 

.04 

.063 

.063 

-.187 

-.187 

■ 

.309 

-.426 

-.426 

.536 

.536 

.08 

.125 

.125 

-.368 

-.368 

1 

.588 

-.771 

-.771 

.905 

.905 

.12 

.187 

.187 

-.536 

-.536 

1 

.809 

-.969 

-.969 

.992 

.992 

.16 

a 249 

.249 

-.685 

-.685 


.951 

-.982 

-.982 

.771 

.771 

.20 

.309 

.309 

-.809 

-.809 

1.000 

1.000 

-.809 

-.809 

.309 

.309 

.24 

.368 

.368 

-.905 

-.905 

.951 

.951 

-.482 

-.482 

-.249 

-.249 

.28 

.426 

.426 

-.969 

-.969 

.809 

.809 

-.063 

-.062 

-.729 

-.729 

.32 

*82 

.482 

-.998 

-.998 

.588 

.588 

.368 

.368 

-.982 

-.982 

.36 

..36 

.536 

-.992 

-.992 

.309 

.309 

.729 

.729 

-.930 

-.930 

.40 

.588 

.588 

-.951 

-.951 

.000 

.000 

.951 

.951 

-.588 

-.588 

.44 

.637 

.637 

-.876 

-.876 

-.309 

-.309 

.992 

.992 

-.063 

-.063 

.48 

.685 

.685 

-.771 

-.771 

-.588 

-.588 

.844 

.844 

.482 

.482 

.52 

.729 

.729 

-.637 

-.637 

-.809 

-.809 

.536 

.536 

.876 

.876| 

.56 

.771 

.771 

-.482 

-.482 

-.951 

-.951 

.125 

.125 

.998 

.998| 

.60 

.809 

.809 

-.309 

-.309 

-1.000 

-1.000 

-.309 

-.309 

.809 

.809: 

.64 

.844 

.844 

-.125 

-.125 

-.951 

-.951 

-.685 

-.685 

.368 | 

.368! 

.68 

.876 

876 

.063 

.063 

-.309 

-.809 

-.930 

-.930 

187 

-. 187 | 

72 

.905 

.905 

.249 

.249 

-.588 

-.588 

-.998 

-.998 

-.685 

-.685 j 

.76 

.930 

.930 

.426 

.426 

-.309 

-.309 

-.876 

-.876 

-.969 

-.969 : 

.80 

.951 

.951 

.588 

.588 

.000 

.000 

-.588 

-.588 

-.951 

-.951 | 

.84 

.969 

.969 

.729 

.729 

' .309 

.309 

-.187 

187 

-.637 

-.637 

.88 

.982 

.982 

.844 

.844 

.588 

.568 

.249 

.249 

-.125 

-.125 

.92 

.992 

.992 

.930 

.930 

.809 

.809 

.637 

.637 

.426 

.426 

.96 

.998 

.998 

.982 

.982 

.951 

.951 

.905 

.905 

.844 

.844 

1.00 

1.000 

1.000 

1.000 

1.000 

1.000 

1 . ^0 

1.000 

1.000 

1.000 

1.000 
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TABLE XV.- COMPARISONS OF COMPUTED UNCOUPLED BENDING MID TORSION NATURAL 


FREQUENCIES WITH REFERENCE SOLUTIONS FOR A CLAMPED, NONROTATING, 
UNIFORM BEAM WITH A TIP NEIGHT* (CASE 2) 


Mode 

Reference frequency,* 5 
cpa 

Computed frequency, 
cpa 

Vertical 

266 

270 

Vertical 

2 837 

2 817 

Torsion 

5 506 

5 505 

Vertical 

8 847 

8 823 

Torsion 

16 517 

16 493 


•The cip weight, as represented, affects only the bending natural 


frequencies. 

^See reference 10. 


TABLE XVI.- COMPARISONS OF COMPUTED UNCOUPLED BENDING AND TORSION 
NATURAL FREQUHICIES WITH EXACT SOLUTIONS POP. A NONROT.' TING, 
SIMPLY SUPPORTED, UNIFORM BEAM (CASE 3) 


Mode 

Computed frequency, 
cptn 

Exact frequency, 
cpm 

Vertical 

1 711 

1 711 

Vertical 

6 846 

6 846 

To* sion 

11 004 

11 012 

Vertical 

IS 402 

15 403 

Torsion 

21 966 

22 023 

Vertical 

2? 380 

! 27 382 

Torsion 

32 841 

33 035 

Vertical 

42 779 

42 785 

Tor s ion 

43 582 

44 047 

Torsion 

54 158 

55 058 
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Figure 2.- Sign conventions for center-of-gravity and 
shear center offsets from pitch axis. 


L^i 



(a) Vertical plane. 



(b) Horizontal plane. 


X 

(c) Torsion plane. 



Figure 3.- Sign conventions for deflections, slopes, shears, and moments 
associated with a blade segment. Arrows indicate positive directions. 
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(b) Forces and aments dependent on deformation. 
Figure 4.- Concluded. 










x(<4 > 




(b) Illustration of natural frequency iteration technique. 

Figure 6.- Determination of natural frequencies using the auxiliary function. 
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Figure 7.- 


Computer program flew chart and logic. 
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Figure 8.- Illustration of program optional output (mode shape plot-*. 
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Figure 9.- Illustration of program optional output (plot of variation 
of natural frequency with rotor speed). 
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Experiment (ref. 12) 
Blade pitch at 0.75R, deg 
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(a) First bending natural frequency. 


Analysis 

Experiment (ref. 12) 



(b) Second bending natural frequency. 






Experiment (ref. 12) 
Blade pitch angle at 0.75R, deg 



(c) Fit st torsion natural frequency. 
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